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Abstract 


Radio  waves  that  pass  through  the  earth's  ionosphere  travel  more  slowly  than 
their  free  space  velocity  due  to  the  group  path  delay  of  the  ionosphere.  This  group 
path  delay,  directly  proportional  to  the  total  electron  content  of  the  ionosphere,  can 
be  an  important  source  of  error  to  VHF,  UHF,  and  L-band  satellite  detection  radars 
and  satellite  navigation  systems.  In  this  report,  the  current  state  of  knowledge  of 
ionospheric  total  electron  content  is  outlined,  with  special  emphasis  placed  on  the 
North  Atlantic  region  of  the  world  due  to  NATO  special  requirements  in  this  region. 

A  numerical  model  of  total  electron  content,  valid  over  the  European  continent  under 
certain  conditions,  is  presented  for  systems  engineering  use  for  an  average  back¬ 
ground  total  electron  content  correction.  Typical  values  of  total  electron  content 
are  also  given  at  various  locations  in  the  high,  middle,  and  equatorial  latitudes.  If 
the  results  presented  here  seem  incomplete,  it  is  only  because  the  state  of  knowledge 
of  the  total  electron  content  parameter  is  still  incomplete.  With  more  observational 
data  being  taken  at  many  locations,  an  over -all  satisfactory  picture  of  the  world-wide 
behavior  of  this  important  parameter  is  beginning  to  emerge. 
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TOTAL  ELECTRON  CONTENT 
STUDIES  OF  THE  IONOSPHERE 

1.  Introduction  -  Importance  of 
the  Total  Electron  Content  Parameter 

J.A.  Klobuchar  and  J.  Aarons 

M.  INTRODUCTION 

The  total  electron  content  (TEC)  of  the  earth's  ionosphere  is  defined  as  the 
total  number  of  free,  thermal  electrons  in  a  unit  area  column  of  ionosphere  from 
the  ground  to  a  height  well  above  the  peak  of  ionization,  to  at  least  1000  km.  This 
TEC  is  usually  described  as  a  vertical  column  one  square  meter  in  area.  Typical 
TEC  values  range  from  101®  to  10^  el/m^,  but  the  exact  value  is  a  function  of 
many  variables.  Some  of  the  variables  which  influence  TEC  are  geographic  location, 
local  time,  season,  solar  EUV  flux,  and  magnetic  activity.  In  a  given  day  at  one 
location,  changes  in  TEC  of  10  to  1,  or  greater,  are  typical  in  some  seasons. 

Why  is  the  TEC  parameter  studied?  What  importance  does  it  have?  First, 
experimentally  it  has  been  a  relatively  simple  parameter  to  measure.  The  earliest 
TEC  results,  obtained  in  the  late  1950s,  showed  the  first  temporal  and  seasonal 
behavior  of  the  topside  of  the  ionosphere.  Comparisons  were  later  made  of  the 
ratio  of  the  topside  to  the  bottomside  content.  The  equivalent  slab  thickness,  the 
ratio  of  TEC  divided  by  the  density  at  the  peak  of  the  F-region,  also  was  studied. 

Slab  thickness  is  the  equivalent  thickness  the  ionosphere  would  have  if  it  had  a 
constant  density  equal  to  the  peak  density  throughout  its  entire  height.  The  slab 
thickness  generally  ranges  from  100  to  400  km  in  height. 
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The  chief  scientific  value  of  TEC  today,  lies  in  its  simplicity  of  measurement, 
providing  a  suitable  satellite  source  is  available.  TEC  provides  continuity  of 
F-region  measurements  through  magnetically  disturbed  periods.  TEC  measure¬ 
ments  have  the  potential  use  as  another  routine  parameter  describing  the  F2  region 
at  many  locations. 

What  engineering  importance  does  TEC  have?  Any  operating  or  potential  sys¬ 
tem  which  involves  radio  waves  propagating  through  the  ionosphere  and  which 
requires  time  delay  measurements  to  an  accuracy  of  the  same  order  as  the  iono¬ 
sphere  time  delay  errors  requires  a  knowledge  of  TEC  to  correct  for  these  errors. 
The  ionosphere  produces  a  retardation  in  the  velocity  of  the  information  carried  on 
a  radio  wave  called  group  path  delay.  This  group 'path  delay  produces  timing  errors 
in  radar  systems  and  navigation  systems  which  must  traverse  the  ionosphere.  The 
group  path  delay  timing  error  is  directly  proportional  to  the  TEC. 

The  group  path  delay  timing  error  can  be  important  depending  upon  the  fre¬ 
quency  the  system  uses-and  the  amount  of  the  ionospheric  TEC.  For  VI1F  and  higher 
the  relationship  bet  ween  time  delay  dt  and  TEC  is: 

dt  =  -^r  TEC  (sec) 
c  f2 

where 

c  =  velocity  of  light  =  3  x  108  mr/sec 

f  =  system  operating  frequency  in  Hz 

TEC  =  total  electron  content  in  el/m^ 

A  typical  time  delay  for  a  radar  operating  at  430  MHz  might  range  from  10  to  1000 
nsec  for  a  target  at  the  zenith.  Figure  1-1  shows  a  plot  of  time  delay  versus  fre¬ 
quency  for  several  TEC  values.  At  25  deg  elevation  these  time  delays  would-be 
approximately  doubled  due  to  the  greater  TEC  encountered  through  an  oblique  iono¬ 
spheric  path.  Figure  1-2  shows  the  factor  by  which  vertical  TEC  values  must  be 
multiplied  to  obtain  TEC  time  delay  values  at  other  elevation  angles.  Typical  mean 
ionospheric  heights  range  from  300  to  400  km  for  this  conversion  ratio.  This  does 
not,  of  course,  correct  for  temporal  or  geographic  gradients  of  TEC  which  produce 
an  additional  difference  between  vertical  and  oblique  TEC  values. 

In  order  to  determine  the  time  delay  correction  for  an  object  at  amarbitrary 
height,  geographic  location  and  time  in  the  ionosphere,  a  complete  description  of 
the  electron  density  height  profile  for  all  times  and- locations  must  be  known.  This 
task  is  well  beyond  the  scope  of  this  report.  Here,  only  the  status  of  measure¬ 
ments  of  TEC,  the  total  vertically  integrated  electron  density  up  to  a  height  of  at 
least  1000  km,  will  be  discussed.  First,  some  comments  will  be  made  on  the 
existing  TEC  data  base,  along  with  prospects  for  improving  this  data  base.  Next, 
the  approximations  which  are  a  part  of  the  TEC  data  base  will  be  discussed.  These 
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ELEVATION  ANGLE  AT  SUAfACE  (EL) 


Figure  1-2.  Multiplication  Factor  to 
Convert  Vertical  Electron  Content  to 
Slant  Content  for  Various  Choices  of 
Mean  Ionospheric  Height 


Figure  1-1.  Time.Delay  Due  to  the  Iono¬ 
sphere  Versus  Frequency  for  Several 
Values  of  Total  Electron  Content 


sections  will  be  followed  by  a  brief 
description  of  efforts  to  date  in  TEC 
modelling.  Finally,  in  the  data  section, 
tables  and  graphs  of  TEC  averages  for 
several  locations  are  presented  along 

with  a  discussion  of  TEC  data  obtained  from  the  three  world  ionospheric  regions- 
the  equatorial,  midlatitude,  and  polar  regions.  A  brief  summary  dif-  usses  possible 
improvements  in  TEC  prediction  capability  in  the  near  term,  ulU  -ate  limitations 
to  TEC  data  and  possible  experimental  work  which  will  improve  TEC  prediction 
accuracy. 


1-2.  SOME  COMMENTS  ON  THE  EXISTING  DATA  BASE 

Measurements  ofothe  total  electron  content  of  the  ionosphere  have  been  made 
since  the  late  1950s.  The  earliest.measurements  were  made  by  monitoring  the 
polarization  twist  of  vhf  lunar  reflected  radio  waves.  Following  the  advent  of 
artificial  earth  satellites  the  Faraday  polarization  rotation  and  differential  Doppler 
effects  were  used  to  determine  TEC  from  VHF  radio  waves  transmitted  from  certain 
satellites.  The  differential  Doppler  and  differential  Faraday  effects  have  also  been 
used  to  determine  the  height  profile  of  electron  density  by  means  of  rocket  probes. 

Numerous  ionospheric  observatories  participate  in  taking  TEC  data  from  one 
or  more  artificial  earth  satellites  but  little  has  been  done  to  integrate  the  results 
from  the  various  independent  observers  and-to  attempt  to  produce  a  world-wide 
picture  of  ionospheric  total  electron  content  from  the  measurements  taken  from  these 
satellites.  In  fact,  it  may  not  be  possible  to  compile  an  over-all  picture  of  TEC 
from  these  data  as  the  available  data  from  many  of  these  stations  consists  of  the 
results  of  a  few  satellite  passes  a  day  averaged  over  a  few  month  period.  These 
seasonal  averages  were  compiled  using  different  times  and  lengths  of  seasons  for 
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the  various  observers.  As  a  result,  combining  these  seasonal  averages  from  the 
various  stations  encounters  the  difficulty  that  even  the  mean  seasonal  picture  may 
suffer  from  these  sampling  differences.  No  day  to  day  variability  can  be  determined 
from  the  mean  seasonal  data. 

One  exception  to  the  rather  bleak  picture  of  the  state  of  the  available  TEC  data 
obtained  from  low  orbit  satellites  is  in  the  European  continent  where  Amayenc, 

Bertin  and  Papet-Lepine  (1971)  combined  TEC  data  from  near  Paris,  France,  with 
that  of  Liszka's  (1967)  data  taken  at  Kiruna,  Sweden,  to  produce  mean  seasonal 
maps  of  TEC  various  local  mean  time  from  39  to  72  deg  North  latitude.  These  TEC 
contours,  graciously  supplied  by  F.  Bertin  for  inclusion  in  this  report,  are  included 
in  the  data  section. 

Geostationary  satellite  polarization  data  are  much  better  with  regard  to  data 
time  sampling,  but,  unfortunately,  the  geographic  distribution  of  the  available  data 
is  severely  restricted  by  the  lack  of  geostationary  satellites  having  suitable  VHF 
transmissions.  Most  of  the  available  total  electron  content  data  from  geostationary 
satellites  have  been  taken  from  the  North  American  continent  and  from  the  Australia- 
New  Zealand  area,  with  data  also  available  from  Japan,  Hawaii,  and  Wales.  The 
European  continent,  until  late  1971,  has  not  had  a  suitable  geostationary  satellite 
available  for  long  periods  of  time  to  make  ionospheric  measurements. 

Perhaps  the  main  point  concerning  availability  of  TEC  data  is  that  there  is  not 
now  -  nor  is  there  likely  to  be  in  the  near  future  —  sufficient  T-EC  data  from  a  large 
enough  number  of  stations  to  construct  a  world-wide  TEC  rr,  '/del  directly  from  TEC 
data  alone.  Therefore,  any  TEC  model,  except  ore  which  covers  only  a  limited 
geographic  region,  must  lean  heavily  on  other  existing  models  of  the  F2  region,  A 
world-wide  foF2  prediction  model  was  made  using  data  from  more  than  50  ionosondes 
for  the  sunspot  minimum  conditions  and  from  over  100  ipnosondes  during  the 
International  Year  sunspot  maximum.  Direct  TEC  measurements  have  little  hope  of 
being  available  from  this  number  of  stations.  Therefore,  models  of  TEC  must 
necessarily  rely  heavily  oa  existing  foF2  data  maps. 
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2.  TEC  From  the  Faraday  Effect 


J.A.  Klobuchar  and  M.  Mendillo 


2-1.  INTRODUCTION 


Measurements  of  the  ionospheric  total  electron  content  (TEC)  are  being  made 
at  several  mid-latitude  sites  by  monitoring  the  Faraday  rotation  of  VHF  radio  waves 
from  geostationary  satellites.  The  correct  conversion  of  the  polarization  data  to. 
equivalent  vertical  TEC  values  is  mainly  dependent  upon  the  choice  of  the  M factor 
which  appears  in  the  Faraday  equation.  A  recent  paper  by  Smitli  (1970)  concludes 
that  large  errors  may  result  if  a  constant  W  factor  is  used  throughout  the  course 
of  a  day.  Here  we  present  the  results  of  calculations  made  of  the  total  Faraday 
rotation  and  the  "true"  number  of  electrons  along  a  path  from  a  mid-latitude  station 
to  a  geostationary  satellite.  Typical  summer  and  winter,  day  and  night  conditions 
were  used  to  determine  how  the  H  factor,  the  polarization  twist  (£)  and  the  TEC 
change  as  a  function  of  height. 

2-2.  THE  FARADAY  METHOD  OF  TEC  MEASUREMENTS 

The  amount  of  total  polarization  twist,  Q ,  is  related  to  the  total  electron  content 
by  the  well  known  relation: 


Q  =  -jp-  J B  cos  0  sec  X  N  dh  (radians) 


(2-1) 
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where: 

K  =  2.36  x  10' 
f  =  frequency  in  Hertz 
B  =  magnetic  field  strength  in  gammas 
9  =  propagation  angle 
X  =  zenith  angle 

N  -  local  electron  density  in  el/M2. 

In  practice,  the  measured  amount  of  polarization  twist,  Q,  is  converted  to  an 
equivalent  total  vertical  electron  content  by  removing  B  cos  i?  sec  X  from  under  the 
integral  sign  and  replacing  it  with  a  mean  value.  Then: 


Q  = 


7w/n 


dh 


(2-2) 


where  "B"cos-5~secAr*  =  H  is  computed  in  the  following  manner.  A  typical  N(h) 
profile  is  assumed  and  calculations  of  the  mean  value  Ivf  are  generally  found  by 
equating: 


M 


.J 


B  cos  6  secX  N  dh 


/  N 


(2-3) 


dh 


The  crux  of  this  discussion  lies  in  the  limits  of  the  integration  in  both  the  numerator 
and  the  denominator  in  Eq.  (2-3).  At  least  one  author,  Smith  (1970),  has  pointed 
out  that  there  is  a  large  diurnal  change  in  Ivf  when  polarization  observations  are 
made  from  signals  transmitted  from  a  geostationary  satellite.  He  used  N(h)  from 
the  Thomson  scatter  facility  at  Arecibo  to  heights  of  wlOOO  km,  the  Angerami  and 
Thomas  (1964)  exospheric  electron  density  model  above  that  height,  and  integration 
limits  from  the  ground  to  the  satellites  in  order  to  determine  the  diurnal  changes  in 
M. 

In  the  calculations  which  follow  we  will  show  that  an  M  derived  from  Eq.  (2-3) 
is  in  error  if  the  integration  is  carried  out  to  the  satellite  height.  The  integration 
of  both  the  numerator  and  the  denominator  in  Eq.  (2-3)  should  only  be  carried  out 
to  heights  where  S3  is  no  longer  measurable,  If  this  point  is  ignored,  the  resultant 

J  N  dh  will  consist  of  both  the  TEC  actually  responsible  for  the  measured  amount  of 

rotation  (£),  plus  an  additional  Jn  dh  along  the  ray  path  where  a  nonmeasurable 
amount  of  rotation  occurs.  In  summary,  then,  Eq.  (2-3)  should  be  rewritten, 
using  Eq,  (2-1)  to  read: 

,2 

(2-4) 


TX  = 


r 


Q 


K 


dh 


7 


where  Q  is  the  measured  amount  of  rotation,  and  the  upper  limit  to  the  integral 
corresponds  to  the  height  above  which  the  remaining  amount  of  rotation  is  less  than 
the  absolute  experimental  error. 


2-3.  ABSOLUTE  ACCURACY  OF  MEASUREMENTS  OF  POLARIZATION  ROTATION 

Before  the  results  of  the  calculations  of  Q  and  TEC  are  presented  a  discussion 
of  the  absolute  measurement  accuracy  of  the  experimental  quantity  in  the  Faraday 
technique,  (Q)  is  in  order.  Measurements  of  the  polarization  twist  of  VHF  radio 
waves  from  geostationary  satellites  are  usually  carried  out  by  means  of  rotating  or 
spinning  yagis.  The  position  of  the  null  in  the  arriving  signal  plane  is  compared 
either  mechanically  or  electrically  with  some  reference  polarization;  for  example, 
the  local  vertical.  Some  workers  use  a  pair  of  right  hand  and  left  hand  circularly 
polarized  helices,  with  a  dual  phase  locked  receiver,  and  compare  the  phase 
difference  at  the  two  receiver  outputs.  The  simple  yagl  or  helical  antennas 
generally  used  to  measure  the  relative  polarization  with  respect  to  a  local  reference 
can  determine  this  angle  on  a  relative  basis  to  within  approximately  +  5  deg,  if 
careful  calibration  of  the  receiving  system  is  carried  out  and  maintained. 

A  correction  also  must  be  made  for  the  angle  that  the  transmitted  signal  from 
the  satellite  would  make  with  respect  to  the  local  vertical  in  the  absence  of  an  iono¬ 
sphere.  The  orientation  that  the  transmitted  VHF  wave  makes  with  respect  to  the 
satellite  spin  axis,  or  to  some  known  reference,  is  usually  not  known  precisely, 
since  measurements  of  this  transmitted  orientation  are  not  made  before  launch. 
Consequently,  resort  must  be  made  to  calculations  of  this  angle  based  upon  the 
satellite  VHF  antenna  and  feed  characteristics.  In  the  case  of  ATS- 3,  an  independent 
determination  of  this:angle  was  made  by  measuring  the  total  twist  of  lunar  reflected 
signals  at  a  high  VHF  frequency  at  a  time  when  the  moon  was  nearly  directly  behind 
the  ATS-3  satellite,  Klobuchar  (1969).  By  knowing  the  absolute  amount  of  twist  on 
the  lunar  reflected  signal  and  relative  amount  of  measured  twist  along  the  ATS-3 
path,  the  transmitted  polarization  of  the  VHF  signal  on  the  ATS-3  was  determined. 

The  accuracy  of  either  the  calculated  or  the  experimental  methods  of  determining 
this  transmitted  polarization  is  certainly  no  better  than  +5  deg.  In  addition  the 
satellite  spin  axis  may  have  some  diurnal  precession  of  a  few  degrees  about  the 
nominal  polar  axis.  In  conclusion,  we  estimate  the  total  polarization  twist  measure¬ 
ment  error  due  to  all  sources,  at  least  for  the  modest  antennas  used  by  most  workers, 
to  be  not  better  than  +10  deg. 


24.  MODEL  CALCULATIONS- ASSUMPTIONS  AND  DATA  USED 


In  order  to  see  how  Q  and  J  N  dh  vary  as  a  function  of  height  along  the  path 
toward  a  geostationary  satellite,  we  carried  out  integration  of  12  and  N  from  a  mid¬ 
latitude  station  at  Hamilton,  Massachusetts,  (42.6°N,  70. 8°W)  to  a  geostationary 
satellite  on  the  station's  meridian.  We  used  N(h)  data  up  to  a  vertical  height  of 
1000  km,  obtained  from  the  Thomson  scatter  radar  at  Millstone  Hill,  kindly  provided 
by  J.V,  Evans  of  the  MIT  Lincoln  Laboratory.  N(h),  above  a  vertical  height  of 
1000  km,  was  derived  from  the  Angerami  and  Thomas  (1964),  (A&T)  diffusive 
equilibrium  exospheric  model.  The  calculations  of  polarization  twist  and  TEC  as  a 
function  of  height  along  the  straight  line  path  from  the  satellite  to  the  station  were 
made  for  four  conditions,  summer  and  winter  day  and  night. 

To  calculate  the  electron  density  at  some  point  along  a  field  line  the  A&T 
exospheric  model  requires  a  temperature,  composition,  and  electron  density  at  the 
base  of  the  exosphere  along  the  same  field  line,  here-taken  as  1000  km.  In  order 
to  compute  the  electron  density  along  the  straight  line  path  to  the  satellite  these 

parameters  are  required  for  the  A&T 
model  over  a  substantial  latitude  range. 
Figure  2-1  shows  that  the  minimum  L 
reached  by  a  ray  from  Hamilton, 
Massachusetts,  to  the  geostationary 
satellite  on  the  station  meridian  is 
L  =  2.3,  which  corresponds  to  an  in¬ 
variant  latitude  at  1000  km  (A1)  of  44  deg. 
L  =  4.6  (A*  =60  deg)  is  taken  as  the 
plasmapause  outside  of  which  we  take  a 
constant  electron  ^density  of  200  el/cm3. 
The  electron  density  at  the  1000  km  base 
level  of  the  exosphere  at  the  latitude  of 
Millstone  Hill  (A1  =  55  deg)  is  taken  from 
the  Thomson  scatter  profiles  provided  by 
Evans.  The  electron  density  over  the 
remainder  of  the  44  to  60  deg  latitude  range  is  taken  from  the  gradients  given  by 
Brace,  Mayr  and  Reddy  (1967)  fitted  to  the  Millstone  Hill  values  at  55  deg. 

Electron  temperatures  over  the  same  latitude  range  were  obtained  by  using  the 
temperature  gradients  of  Brace  et  al,  fitted  to  the  actual  Te  values  measured  at 
1000  km  at  Millstone  Hill,  also  "kindly  provided  by  Evans.  In  the  A&T  model  of 
diffusion  along  a  field  line  the  ion  and  electron  temperatures  are  assumed  to  be 
equal  to  the  1000  km  value.  Ion  composition  was  chosen  on  the  basis  of  measured 


Figure  2-1.  Path  of  the  Straight  Line 
Ray  From  a  Geostationary  Satellite  to 
Sagamore  Hill,  Hamilton,  Massachusetts 
(42.6°N.,  70. 8°W. ),  The  satellite  is  on 
the  station  meridian 
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values  at  Arecibo  (A1  =  33  deg),  Prasad  (1970),  and  estimated  values  at  55  deg 
fitted  with  a  linear  change  between  these  two  latitudes. 

2-5.  MODEL  CALCULATIONSEXOSPHERIC  CONTRIBUTION  TO  THE  TEC 

Results  of  the  calculations  of  Q  and  TEC  from  the  satellite  to  the  station  for 
the  four  cases  being  considred  are  shown  in  Figure  2-2.  Both  Q  and  TEC  above 
the  ordinate  height  are  plotted  as  a  function  of  height.  The  values  above  1000  km 

in  vertical  height,  shown  in  Figure  2-2, 
are  from  the  A&T  model  calculations. 
-  -  Below  1000  km  Q  and  TEC  are  from 

*«T«a*n  •iifiiMMT 

«>'  — —  Evans'  Ng(h)  profiles.  All  values  were 


TV  computed  along  the  slant  path  at  41  deg 

\  \  t  \  \ 

\  I  «N\  elevation  angle,  though  the  ordinate  is 

the  more  commonly  used  vertical  height. 

i  t . .  l\ ^  J  .  .4.^ . -II  ■■■.„ ;  The  point  where  10  deg  of  rotation  occurs, 

l  ",mn  \  *•""  1  which  is  within  the  estimated,  error  of 

c  J  “■”**"  measurement  accuracy  of  most  observers, 

inm.  is  indicated  on  the  abscissa  by  a  small 

\  \  \  \  arrow.  Note  that  above  this  point,  where 

"**:  \V  only  10  deg  of  rotation  remains,  there  is 

y,  still  a  significant  TEC.  Table  2-1 

5*  f.  indicates  both  the  amount  and  the  percent 

TEC  remaining  along  the  path  to  the 

Figure  2-2.  Polarization  Twist  (Q,  in  satellite  above  the  10  deg  point,  and  the 

Radians)  and  Total  Electron  Content  men 

(TEC.  in  el/m2)  Along  the  Slant  Path  TEC  above  1000  ^  vertical-height. 

From  a  Geostationary  Satellite  to  a  Mid-  As  Table  2-1  shows,  there  is  still  a 

Latitude  Station  .  ...  .  .  ,  ,  .. 

significant  percentage  of  TEC  above  the 

point  where  the  remaining  Faraday  rota¬ 
tion  is  within  the  measurement  error.  Hence,  the  Faraday  rotation  effect  does  not 
respond  measurably  to  this  portion  of  the  "true"  total  electron  content. 

Table  2-2  gives  the  results  of  computing  M  according  to  Eq.  (2-4),  including 
bothi?  and  TEC  up  to  three  heights;  1000  km,  the  point  above  which  the  remaining 
rotation  is  less  than  10  deg,  and,  the  satellite  height.  Also  shown  in  Table  2-2 
are  the  percentage  changes  for  each  condition  from  a  mean  of  the  four  cases. 

The  results  in  Table  2-2  show  that  the  day  to  night  change  in  the  M  factor  is 
smallest  for  "KT  derived  by  using  Q  and  TEC  up  to  1000  km.  If  "M  is  obtained  by 
choosing  the  integration  limits  in  Eq.  (2-4)  at  the  height  of  10  deg  of  remaining 
rotation  the  day  to  night  change  is  still  rather  small.  The  worst  case  is  in  winter. 


—  to*  io"* 

“  «*ur>  t  ii*n  i*  »v  «»uri  .it  on  i* 
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Figure  2-2.  Polarization  Twist  ( Q ,  in 
Radians)  and  Total  Electron  Content 
(TEC,  in  el/m2)  Along  the  Slant  Path 
From  a  Geostationary  Satellite  to  a  Mid- 
Latitude  Station 
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Table  2-1.  Amount  and  Percent  TEC  Remaining  Along  Path  to  Satellite  Above 
10  deg  Point  and  TEC  and  S2  Above  1000  km  Vertical  Height 


TEC 
above 
1000  km  . 

%  TEC 
above 
1000  km 

a 

above 
1000  km 

%Q 
above 
1000  km 

TEC 

above 

10  deg  point 

%  TEC 
above 

10  deg  point 

Winter  Day 

5.19 

12.2 

O 

o 

3.4 

3.72 

8.8 

Winter  Night 

4.53 

30.8 

47° 

13.7 

2.73 

18.6 

Summer  Day 

5.35 

18.6 

51° 

6.9 

3.56 

12.4 

Summer  Night 

2.91 

31.3 

25° 

11.6 

2.34 

25. 4 

Table  2-2.  Results  of  Computing  M  According  to  Eq.  (2-4) 


M1000 

M10° 

^total 

^1000 

(Pert 

M10deg  |  ^total 
:ent  Changes). 

Winter  Day 

60236 

59251 

52906 

2.9 

4.6 

13.9 

Winter  Night 

57062 

53434 

42508 

-2.5 

-5.6 

-8.5 

Summer  Day 

58353 

56934 

48577 

CO 

© 

1 

-0.5 

4.6 

Summer  Night 

58475 

56908 

41854 

0 

0.5 

-9.9 

where  a  10  percent  day  to  night  error  is  calculated.  However,  if  "M  is  calculated 
using  integration  limits  to  the  height  of  the  satellite,  the  day  to  night  "M  change  is 
22.4  deg  in  winter  and  14.5  deg  in  summer.  The  meanM  is  18  percent  less  than 
the  M  calculated  to  the  10  deg  point  and  20.6  percent  less  than  M  computed  using 
integration  limits  of  1000  km.  The  difference  between  the  mean  M  for  1000  km 
and  for  the  1 0  deg  height  point,  however,  is  only  3.  2  percent. 


2  6.  DISCUSSION 

We  have  shown  examples  of  four  cases  of  calculations  of  iJand  TEC  calculated 
up  to  a  geostationary  satellite  from  a  mid-latitude  station.  The  four  cases  were 
chosen  to  represent  ionospheric  diurnal  maximum  and  minimum  conditions  where 
the  contribution  from  the  exosphere  might  be  expected  to  be  a  minimum  and  a  maxi¬ 
mum  of  the  total  TEC  respectively.  The  A&T  model  used  is  a  constant  temperature 
model  along  field  lines  and  requires  a  knowledge  of  temperature,  composition  and 
electron  density  at  the  base  level  over  a  44  to  60  deg  latitude  range.  Because  of  the 
several  parameters  of  the  model,  which  itself  is  an  approximation,  these  calcula¬ 
tions  may  not  reflect  the  actual  exospheric  densities.  They  do  show,  however,  that 
the  polarization  twist  predominantly  occurs  in  the  ionosphere,  with  only  from  3  to 
14  percent  occurring  above  1000  km. 
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In  summary,  the  Faraday  technique,  when  used  with  VHF  signals  transmitted 
from  geostationary  satellites,  remains  a  good  measure  of  the  ionospheric  total 
electron  content;  but  it  is  not- a  satisfactory. measure  of  the  exospheric  electron 
content.  Cnanging  the  M  factor  in  an  attempt  to  account  for  the  electrons  above  the 
point  where  they  are  measured  by  the  Faraday  effect  is  misleading,  we  believe,  and 
can  be  in  serious  error.  We  maintain  that  a  more-satisfactory  method  of  presenting 
TEC  results  obtained  fromVHFsignals  from  geostationary  satellites  is  to  use  an  M 
factor  calculated  by  using  typical  profiles  up  to  a  standard  height  of  1000  km  and  to 
include  any  exospheric  TEC  contribution  as  an  additive  constant. 


2-7.  DETERMINING  DIURNAL  IN  MEAN  FIELD  HEIGHT 

After  this  section  of  the  report  was  written,  a,paper  by  J.  E.  Titheridge  (1972) 
was  published  in  which  he  used  a  wide  range  of  model  ionospheres  and  satellite  and 
station  locations  to  determine  diurnal  changes  in  the  mean  field  height.  He  con¬ 
cluded  that  a  fixed  mean  height  of  420  km  gives  a  resultant  accuracy  of  plus  and 
minus  5  percent  in  TEC  up  to  a  mean  height  of  approximately  2000  km  under  most 
conditions.  Our  work  is  in  agreement  with  his  conclusions,  which  are  more  general 
in  geographic  "coverage. 

Another  point  concerning  the  conversion  of  the  Faraday  polarization  twist  to 
TEC  has  been  raised  by  R.  Fritz  (private  communication).  He  pointed  out  that  the 
diurnal  change  in  position  of  a  geostationary  satellite  due  to  a  non-zero  orbital 
inclination  can  produce  a  change  in  the  mean  M  factor,  even  at  the  constant  mean 
height,  of  plus  and  minus  5  percent.  A  knowledge  of  the  precise  orbit  of  the  satellite 
can  be  used  to  correct  for  this. error. 
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3.  Total  Electron  Content  Models 

JA.  Klobuchar  and  R.S.  Allen 


Recently  there  has  been  much. interest  by  systems  engineering  groups  and  by 
other  non-ionospheric  workers  in  the  time  delay  caused  by  the  number  of  free  elec¬ 
trons  in  the  earth's  ionosphere.  These  groups  include  thosewho  require  error 
correction  for  UHF  radars,  proposed  geostationary  satellite  navigation  systems, 
satellite  tracking  stations  using  VHF  beacons,  and  radio  astronomers  who  are 
attempting  to  measure  accurate  positions  and  distances  by  means  of  very  long  base¬ 
line  interferometry.  As  a  result  of  all  this  interest  several  models  of  TEC-time 
delay  have  been  made  by  different  groups.  Most  of  these  models  have  been  con¬ 
structed  using  a  limited  data  base,  hence,  the  applicability  both  in  geographical 
extent  and  in  time  in  the  solar  cycle  is  limited.  A  brief  outline  of  these  models  is 
presented  below. 


3-1  THE  ESSA-NOAA  IONOSPHERIC  PROFILE  MODEL 


The  Environmental  Science  Services  Administration  [Wright,  (1967)],  does 
true  height  analysis  of  ground  based  ionosondes  to  find  the  electron  density  up  to 
the  peak  of  the  F2  region.  Above  the  peak  they  use  the  following  expression: 

N  =  Nmax  •  exp|  {(1  +  g)  •  (1  -  z  -  «'z)} 


where 


=  im 


l  +  g 


<h  •  hmax) 

- IT - 


g  =  0.05 
H  =  scale  height 
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This  resulting  profile  is  accurate  up  to  the  peak  of  the  F2  region  and  is  a  fair  approxi 
mation  to  the  shape  of  the  upper  F  region  above  the  peak.  The  parameter  g,  which 
is  a  measure  of  the  increase  in  scale  height  with  height,  does  have  a  diurnal  and  geo¬ 
graphic  dependence  and  cannot  be  represented  by  a  single  constant.  The  obvious 
disadvantage  to  this  modelcis  the  requirement  for  a  true  height  ionogram  analysis. 
Usually  monthly  median  hourly  ionograms  are  constructed  and  the  true  height 
analysis  is  done  only  once  for  each  hour,  making  a  total  of  24  true  height  analyses 
required  for  each  month. 


3-2.  THE  AFCRL  FIRST  ORDER  MID  LATITUDE  TEC  MODEL 

A  model  of  mid-latitude  slab  thickness  was  made  at  AFCRL  [Klobuchar  and 
Allen  (1970)]  ,  which,  when  combined  with  either  a  predicted  or  an  actual  value  of 
foF2,  gives  a  time  delay.  This  model  is  simply: 

TEC  =  1.24  x  1013  (foF2)2  ^261  +  26  sin  [— +  K  sin  [-ygy—  | 

where 

H  is  the  local  hour  at  the  subionospheric  point  where  TEC  is  desired 

K  =  73  for  local  hours  06  to  19 

K  *  36  for  05  and  20  hours  only 

K  =  0  for  21  to  04  hours 

TEC  is  in  el/M2 

D  is  the  day  of  the  year 

Since  this  model  used  only  northern  mid-latitude  slab  thickness  values  in  its  construe 
tion  its  use  is  limited  to  that  geographic  region. 


3-3  SAMSO  STUDIES 

The  United  States  Air  ForceSpace-and  Missile  Systems  Organization  (SAMSO), 
interested  in  ionospheric  errors  in  a  proposed  L-fcand  satellite  navigation  system, 
funded  three  groups  to  construct  and  test  models  of  the  TEC  over  the  Continental 
United  States  (CONUS),  Puerto  Rico,  Hawaii,  and  Southern  Alaska.  These  studies 
are  now  finished  and  results  of  two  of  the  three  groups  have  been  published  in  final 
form.  The  main  difference  between  the  three  studies  consisted  in  the  approaches  to 
the  modelling  of  the  ionosphere  that  were  taken.  Stanford  University  [iValdman  ana: 
daRosa,  (1971)]  ,  one  of  the  contractors,  constructed  a  TEC  model  based  only  on 
TEC  data.  They  essentially  made  a  polynomial  fit  to  TEC  behavior  as  a  function  of 
solar  activity  and  did  a  Fourier  analysis  of  TEC  diurnal  behavior.  This  was  done 


for  data  from  only  two  stations,  one  was  Stanford,  a  mid-latitude  station  and  the 
other;  Hawaii,  an  equatorial  anomaly  station. 

The  second  SAMSO  contractor,  the  University  of  Illinois  [Rao,  et  al  (1971)] 
used  TEC  data  and  predicted  foF2  values  for  the  ESSA  prediction  program  to 
construct  equivalent  slab  thickness  values  at  one  location,  which  was  then  assumed 
to  be  true  over  a  second  location.  The  ESSA  foF2  prediction  program  then  was  used 
as  a  cornerstone  of  the  data  base  along  with  the  TEC  data  that  were  available  to 
form  predictions  of  TEC  at  another  location. 

The  third  contractor,  the  Applied  Physics^Laboratory  of  the  Johns  Hopkins 
University,  made  an  ionospheric  profile  model  based  upon  ESSA  predicted  foF2 
values  and  theoretical  valuesof  profile  shape  above  the  height  of  the  maximum  of 
the  F2  region.  They  used  TEC  experimental  data  only  in  the  evaluation  of  the  errors 
of  their  model  but  not  in  its  construction. 

34  NASA-DBA  STUDIES 

NASA,  [Bent,  et  al  (1971)]  ,  Goddard  Space  Flight  Center,  recently  had 
completed  for  it  a  study  of  ionospheric  profile  modelling  by  the  DBA  Systems 
Incorporated  Company  of  Melbourne,  Florida.  The  NASA-DBA  model  is  a  profile 
model  constructed  using  the  ESSA  foF2  prediction  program  for  the  peak  density  and 
the  ESSA  MUF  predictions  for  the  height  of  the  ionosphere.  Values  of  scale  height 
above  the  peak  were  determined  from  analysis  of  topside  ionosonde  profiles. 

3-5  COMMENTS  ON  MODEL  STUDIES 

All  of  the  four  F2  region  models  outlined  have  limitations  due  to  the  available 
data  base.  The  most  severely  geographically  limited  model  is  perhaps  the  Stanford 
one,  while  the  one  which  has  the  most  world-wide  applicability,  is  perhaps  the  NASA- 
DBA  model.  Basically,  there  are  at  least  two  different  philosophies  presented  by 
the  four  models  just  completed.  One  approach  was  to  model  TEC  directly  and  only 
to  use  TEC  data  for  this  model.  The  opposite  approach  is  to  model  the  complete 
foF2  profile  with  the  TEC  being  only  a  byproduct  of  the  complete  profile.  These 
different  approaches  meet  different  needs.  The  TEC  only  model  is  of  value  where 
an  adequate  TEC  data  history  exists  to  make  a  direct  model  and  where  only  TEC  is 
required.  That  is,  this  model  is  useful  for  error  correcting  from  targets/satellites 
well  above  most  of'  the  ionization,  say  at  least  1000  km  in  height.  The  profile  model 
is  necessary  if  the  satellite /target  is  within  the  ionosphere  and  the  integrated  elec¬ 
tron  density  only  up  to  the  height  of  the  target  is  required. 


\ 

\ 
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3  6  THE  PENNSYLVANIA  STATE  MARK  I  MODEL 

A  theoretical  model  of  the  ionosphere  has  been  developed  at  the  Pennsylvania 
State  University  [Nisbet  (1971)]  which  attempts  to  describe  the  behavior  of  the 
ionosphere  beginning  with  the  basic  processes  which  control  the  production,  loss 
and  transport  of  electrons.  Experimental  electron  density  values  are  used  only  to 
fit  the  boundary  conditions.  The  model  is  basically  for  the  mid-latitude  ionosphere. 


3-7  LIMITATIONS  OF  PRESENT  MODEL  STUDIES 

The  Stanford  and  Illinois  results  showed  that  the  ionospheric  time  delay  error 
can  be  reduced  by  70-90  percent  over  the  CONUS  with  the  models  they  recently 
completed.  Errors  over.  Hawaii  are  substantially  larger.  No  eastern  longitude  or 
Arctic  latitude  or  southern  hemispheric  data  were  included  in  o.?  tested  against 
their  models.  The  NASA-DBA  profile  model  has  not  been  checked  against  actual 
TEC  data  in  a  statistical  way;  however,  the  monthly  median  predictions  of  TEC 
from  the  NASA-DBA  model  fall  within  the  extreme  Stanford  actual  TEC  monthly 
values  for  the  periods  that  were  tested  in  this  manner. 

All  models  suffer  from  incomplete  data-in  the  polar  regions  and  tne  direct  TEC 
model  suffers  because  of  the  TEC  data  available  which  was  limited  to  the  United 
States. 


3-8  ADVANTAGES  OF  COMBINED  TEC  AND  SLAB  THICKNESS  MODEL 

Using  TEC  obtained  from  polarization  twist  measurements  from  geostationary 
satellites  and  slab  thickness  values  obtained  from  nearby  ionosondings,  continuous 
ionospheric  measurements  including  a  measure  of  the  topside  can  be  made.  Top¬ 
side  satellite-borne  sounders  give  only  limited  temporal  coverage.  The  shape  of 
the  F2  region  is  not  expected  to  vary  as  rapidly  as  the  maximum  density  and  hence 
the  slab  thickness  parameter  should  be  a  good  way  to  combine  foF2  measurements 
or  predictions  with  TEC  measurements  to  make  a  TEC  prediction  at  a  distant  point. 
It  may  also  be  reasonable  to  obtain  a  scale  height  from  slab  thickness  values  and 
make  up  a  profile  model  using  TEC  measurements  and  the  ESSA  foF2  predictions. 
Providing  the  geostationary  satellite  with  a  suitable  vhf  transmitter  is  available  for 
use,  the  TEC  measurements  are  relatively  easy  to  make. 
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4.  Available  TEC  and  SLAB  Thickness  Data 

From  S-66  Satellites 

J.A.  Ktobuchar,  J.P.  Mullen  and  D.R.  Seaman 


4-1  INTRODUCTION 

Many  groups  have  made  TEC  measurements  using  the  Faraday  effect  from  VHF 
beacon  signals  on  the  S-66  (BE-B  and  BE-C)  satellites.  However,  little  effort  has 
been  made  to  use  the  available  data  from  these  stations  to  attempt  to  make  a  model 
of  TEC,  primarily  because  of  the  difficulty  of  combining  data  taken  by  several  workers 
who  used  different  data  collection  and  reduction  techniques,  as  well  as  different 
periods  for  compilation  of  seasonal  average  data.  A  few  studies  of  short  periods  of 
data  have  been  done.  Frihagen  [  JSSG,  (1968)]  made  a  comparison  of  TEC  obtained 
froma  few  satellite  passes  from  several  stations  over  Europe  in  January  1965.  He 
showed  that  similar  latitude  gradients  were  obtained  on  several  satellite  passes  taken 
from  three  different  northern  European  stations.  However,  his  limited  data  set  pre¬ 
cluded  making  meaningful  synoptic  maps  of  TEC.  Houminer  [JSSG,  1968B)]  made 
a  comparison  of  TEC  data  from  three  low-altitude  stations.  The  three  stations  had  a 
similar  diurnal  TEC  variation,  with  the  equatorial  station  having  higher  TEC  values 
but  lower  slab  thickness  values  than  the  other  two  stations.  His  comparison  was  for 
only  one  season.  Klobvchar  and  Aarons  (1968)  combineC'TEC  data  taken  from 
Narssarssuaq,  Greenland,  Hamilton,  Massachusetts,  and  Aricebo,  Puerto  Rico  to 
obtain  the  TEC  and  slab  thickness  latitude  dependence  over  the  5  to  65  deg  North 
latitude  range,  but  only  for  the  March  1966  period.  Galdon  and  Alberca  (1971) 

.studied  the  seasonal  and  solar  cycle  variation  of  mid-day  TEC  data  taken  during 
the  1965  through  1967  period  from  three  European  stations.  They  found  similar 


Preceding  pat* ** 


20 


behavior  with  10. 7  cm  solar  flux  for  a  13  month  running  mean  of  TEC  data  from  all 
three  stations.  After  subtracting  the  13  month  running  mean  solar  cycle  dependence 
from  the  data,  a  large  seasonal  variation  of  TEC  remained  which  changed  with  solar 
flux.  The  magnitude  of  the  seasonal  variation  was  different  for  all  the  three  stations. 
They  used  only  mid-day  TEC  data  to  determine  the  TEC  solar  flux  dependence. 

Several  workers  have  constructed  iso-contours  of  TEC  versus  latitude  and 
local  time  for  one  or  more  seasons.  TEC  values  over  the  range  of  plus  and  minus 
10  to  15  deg  of  latitude  can  be  obtained  in  th'is  manner  from  a  single  station. 

These  contours  are  helpful  in  determining  the  geographic  regions  and  times  where 
certain  anomalies  in  TEC  occur.  All  of  the  available  TEC  contours  are  included 
in  the  data  section  of  this  report. 

Perhaps  the  best  set  of  TEC  contours  was  made  by  Amayenc,  et  al  (1972)  in 
which  they  combined  the  seasonal  contours  of  Liszka  (1967)  with  those  from  their 
own  station  at  Val-Joyeux,  near  Paris,  France  to  produce  TEC  contour  maps  versus 
local  time  over  the  European  sector  from  40  deg  to  72  deg  latitude  for  nine  seasons 
in  the  1964  to  1967  period.  Since  this  geographic  region  is  an  important  one  for 
knowing  the  ionospheric  time  delay  correction  for  a  proposed  VHF navigation  system 
and  for  UHF  satellite  detection  radars,  it  is  fortunate  that  the  work  of  Amayenc  et  al 
and  Liszka  has  been  done  in  this  combined  fashion.  A  separate  section  of  this 
report  details  a  numerical  model  of  TEC,  made  from  their  data,  for  the  year  1965, 
over  the  European  longitude  sector  from  40  to  70  deg  latitude.  This  model  repre¬ 
sents  the  most  complete  direct  TEC  data  set  available  over  Europe,  though  it  is 
only  for  near  sunspot  minimum  conditions.  Several  stations  in  Europe  either  have 
be  i,  or  currently  are  taking  TEC  data,  using  VHF  beacons  on  a  geostationary 
satellite,  but,  thus  far,  little  of  this  data  is  in  finished  form  for  model  studies. 

While  many  papers  in  the  literature  give  details  of  the  behavior  of  TEC  at  a 
particular  latitude,  it  is  probably  not  possible  to  construct  directly-a  complete 
model  of  TEC  from  the  separate  sets  of  published  data.  Fortunately,  the  slab  thick¬ 
ness  parameter  S,  does  not  change  with  latitude  as  much  as  either  TEC  of  Nmax, 
Also,  the  variation  of  slab  thickness  with  sunspot  number  is  considerably  smaller 
than  that  of  either  TEC  or  Nmax.  Thus,  the  slab  thickness  parameter  is  worth 
investigating  for  model  use.  Since  slab  thickness  values  have  been  published  for 
several  stations  having  wide  geographic  distribution,  these  values  have  been  scaled 
off  the  curves  andothe  published  data  are  given  here.  In  a  few  cases  these  values  of 
S  were  made  from  TEC  data  taken  from  observation  of  the  Faraday  effect  on  VHF 
signals  from  a  geostationary  satellite.  A  list  of  tlie  stations  for  which  S  values  are 
available  is  given  in>  Table  4-1  along  with  the  station's  geographic  coordinates,  geo¬ 
magnetic  latitude  and  L  value.  Also,  the  season  for  which  the  data  are  available  is 
indicated  in  the  table.  Figure  4-1  shows  the  behavior  of  S  for  the  winter  and  summ.r 
seasons  for  the  various  stations  for  which  data  were  available.  Note  that  the  daytime 
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values  are  much  greater  in  summer  at  all  latitudes,  than  in  the  winter.,  The  large 
values  of  S  found  in  the  winter  night  at  Honolulu,  particularly  near  sunrise  and  sun¬ 
set,  may  be  due  to  small,  systematic  errors  in  TEC  or  to  relatively  large  horizontal 
gradients  between  the  TEC  measurement  sub -ionospheric  point  and  the  location  of 
the  ionosonde  at  those  times. 

Liszka  (1967)  has  published  iso-contours  of  equivalent  vertical  slab  thickness 
for  northern  Europe  for  three  seasons  during  the  l964  to4965  time  period.  His 
iso-contours  are  given  in  the  data  section  of  this  report.  Also  in  the  data  section 
are  examples  of  comparisons  of  slab  thickness  values  taken  at  widely  separated 
mid-latitude  stations,  indicating  that  the  concept  of  modelling  slab  thickness,  along 
with  a  knowledge  of  foF2,  to  find  TEC,  is  valid.  As  more  TEC  data  become  avail¬ 
able,  particularly  from  monitoring  VHF  radio  waves  from  geostationary  satellites, 
the  behavior  of  slab  thickness  will  become  known  better  than  at  present. 


4-2  STATIONS  WHERE  TEC  DATA  MAY  BE  AVAILABLE 

An  attempt  has  been  made  here  to  list  the  stations  located  throughout  the  world 
where  TEC  data  may  be  available.  These  stations  generally  fall  into  several  cate¬ 
gories  because  of  the  availability  of  appropriate  satellites.  Generally  those  available 
fall  into  two  distinct  sets.  The  earliest  set  of  data  came  from  stations  scattered 
throughout  the  world  who  obtained  TEC  from  various  low-altitude  (approximately 
1000  km)  satellites.  These  satellites  were  the  Transit  series  where  a  few  stations 
obtained  TEC  from  the  differential  Doppler  effect,  and,  later,  the  S-66  satellites 
BE-B  and  BE-C,  where  most  stations  used  the  differential  Faraday  effect  to  obtain 
TEC.  Beginning  in  the  Pacific  sector  in  1964  some  continuous  TEC  data  have  been 
obtained  because  of  the  availability  of  continuous  VHF  transmissions  from  a  geo¬ 
stationary  satellite.  In  1965  the  American  sector  first  had  the  continuous  availability 
of  VHF  signals  from  geostaionary  satellites.  The  European  sector  coverage  has  not 
been  continuous  due  to  the  lack  of  VHF  signals  on  geostationary  satellites  visible  from 
Europe.  Some  European  stations  obtained  continuous  data  during  the  summer  of  1965, 
then  sporadically  b  jginning  again  in  1968.  Thus,  the  available  data  are  limited  in 
various  places,  not  only  because  of  the  obvious  lack  of  stations  in  some  areas  but 
also  due  to  the  absence  of  a  suitable  satellite  of  opportunity  having  VHF  transmissions 
from  which  the  Faraday  rotation  measurement  can  be  made.  TEC  values,  when 
available  from  a  station,  are  much  more  complete  than  those  from  low- orbit  satellites, 
as  they  are  continuous.  The  low  orbit  data  usually  consists  of  a  few  points  per  day 
on  a  combined  seasonal  plot.  Only  in. a  few  instances  are  latitude  gradients  available 
from  the  low  orbit  derived  data. 
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Some  errors  are  inevitable  in  any  list  such  as  this  and  the  author  apologizes  in 
advance  for  any  errors  and  omissions.  They  are  unintentional-.  It  is  simply  not 
possible  to  keep  track  of  the  work  of  all  TEC  experimenters  and  their  data  through¬ 
out  the  world  at  all  times. 

The  stations  have  been  grouped  according  to  three  world  sectors:  the 

American  sector;  2,  the  European-African, sector;  and  3,  the  Asian- Pacific  sector. 
Also  listed  for  each  station  are  the  approximate  geographic;coordlnates  for  each 
station,  the  principal  investigator  or  pereon  otherwise  responsible  for  the  data,  and, 
in  cases  where  it  is  known,  the  availability  of  data  signified  by  S  if  it  consists  of 
TEC  values  derived  from  a  few  passes  per  day  of  a  low  orbit  satellite,  and  G  if  the 
data  came  from  geostationary  satellite  beacon  signals.  The  G  data  are,  of  course, 
continuous  and  presumed  to  be  of  better  quality.  These  stations  are  listed  in 
Table  4-2. 


Table  4-2.  Stations  Where  TEC  Data  May  be  Available 


American  Sector 

Station  Name 

Latitude 

East 

Longitude 

Principal  Investigator 
and  Affiliation 

Data 

Available 

Thule,  Greenland 

77 

291 

J.  Klobuchar,  AFCRL 

G 

College,  Alaska 

65 

212 

G.  Stanley,  U.  of  Alaska 

S,  G 

Baker  Lake,  Canada 
Narssarssuaq.  Greenland 

64 

264 

G.  Swenson,  U.  of  Illinois 

s 

61 

315 

A.  Lundbak,  Danish  Meteorological 

G 

Cold  Bay,  Alaska 

55 

197 

Institute 

K.C.  Yeh,  U.  of  Illinois 

G 

Goose  Bay,  Labrador 

53 

300 

J,  Klobuchar,  AFCRL 

G 

Houghton,  Michigan 

47 

271 

G.  Swenson,  U.  of  Illinois 

S 

Hamilton,  Massachusetts* 

43 

289 

J.  Klobuchar,  AFCRL 

S.W.G 

London,  Ontario 

43 

279 

G.  Lyon,  U.  of  Western  Ontario 

S,  G 

Weston,  Massachusetts 

42 

289 

B.  Relnlsch,  Lowell  Tech. 

S 

University  Park,  Pennsylvania 

41 

282 

Research  Foundation 

W.  Rosa,  Penn.  State  U. 

s,-w 

Boulder,  Colorado 

40 

255 

K,  Davies,  National  Oceanographic 

Ft,  Monmouth,  New  Jersey 

40 

286 

and  Atmospheric  Administration 

H,  Solcher,  Institute  for 

S,  G 

Chesapeake  Beach,  Maryland 

39 

283 

Exploratory  Research 

J.  Goodman,  Naval  Research  Labs 

G 

Blossom  Point,  Maryland 

33 

283 

L.  Blumle,  NASA  Goddard 

S 

Greenbelt,  Maryland 

38 

283 

S.  Rangaswamy,  NASA  Goddard 

G 

Urbana,  Illinois* 

38 

273 

K.C.  Yeh,  University  of  Bllnols 

S,W,G 

Stanford,  California* 

37 

238 

,  A.V.  daRosa,  Stanford  University 

S,  G 

Stanford,  California 

37 

238 

'  V.  Eshleman,  Stanford  University 

China  Lake,  California 

35 

243 

Pioneer 

W,G 

Goldstone,  California 

35 

243 

B,  Mulhall,  Jet  Propulsion 

Laboratory 

G 

Corona,  California 

Los  Angeles,  California 

34 

243 

w. 

-34 

242 

S.  Venkateswaran,  U.  of  California 

G 

College  Station,  Texas 

31 

264 

at  Los  Angeles 

J.  German,  Texas  A  &  M 

s 

Areclbo.  Puerto  Rico 

Kingston,  Jamaica 

18 

293 

A.V.  daRosa,  Stanford  University 

G 

18 

283 

P.  Chin,  U.  of  the  West  Indies 

G 

Mtahuatlan,  Mexico 

Panama  Canal  Zone 

16 

9 

264 

280 

J.  Klobuchar,  AFCRL 

G 

Lima,  Peru 

-11 

283 

A.V,  daRosa,  Stanford  University 

G 

24 


Table  4-2.  (Continued) 


East 

Principal  Investigator 

Data 

Station  Name  .  . 

Latitude 

Longitude 

and  Affiliation 

Available 

Huancayo,  Peru 

-12 

215 

A.  Glesecke,  Geophysical 

S,  G 

Institute  of  Peru 

SanJoae,  Brazil 

-23 

315 

F.  de  Mendonca 

S,W.G<?) 

Tucuman.  Argentina 

-27 

297 

S.  Radies’. la,  7 

S,W 

Port  Stanley 

-51 

303 

S.W 

European  -  African  Sector 

East 

Principal  Investogator 

Data 

Station  Name 

Latitude 

Longitude 

and  Affiliation 

Available 

Tromso,  Norway 

70 

19 

O.  Bratteng  Auroral  Observatory 

S.W 

Kiruna,  Sweden 

65 

21 

L.  Liszka,  Kiruna  Geophysical 

S.W.G 

Observatory 

Kjeller,  Norway 

60 

11 

J.  Frlhagen,  Norwegian  Defense 

S.W 

Research  Establishment 

Rude  Skov.  Denmark 

56 

12 

A.  Lundbak,  Danish  Meteorological 

S 

Institute 

Motcow,  USSR 

55 

37 

J.  Al'pert,  Izmiran 

? 

Kuhlungsbom 

54 

12 

S.W 

Lancaeter.  England 

54 

357 

A.  Hunter,  U.  of  Lancaster 

? 

JodrellBank,  England 

53 

358 

G.  N.  Taylor,  Royal  Radar 

S 

Establishment 

Aberyatwyth,  W»lee» 

52 

356 

L.  Kersley,  U.  of  Wales 

S.  G 

Llndau,  Germany 

52 

10 

•G. .  Hartmann,  -  Max-  Planck 

S.W.G 

Institute  for  A  economy 

Bochum.  Germany 

51 

7 

H,  Kaminski,  Observatory  of  Bochum 

S.W 

Slough.  England 

51 

359 

E.  Golton,  Radio  Research  Station 

Lannlon,  France* 

49 

357 

J.  Papet-Leplne, -National  Center 

S,  G 

for  the  Study  of  Telecommunications 

Breiaach,  Germany 

48 

8 

C.  Munther,  Breisach  Isospheric 

S.W 

Institute 

Florence,  Raly 

44 

11 

P.  Checcacci,  Institute  for  the  Study 

S.W.G 

of  Electromagnetic  Waves  (IROE) 

Val-Joyeux,  France 

44 

4 

F.  Bertln,  National  Center  for 

S.W 

Scientific  Studies  (CNES) 

Graz,  Austria 

41 

15 

R.  Leitlnger,  U.  of  Graz 

S.W 

Tortoaa,  Spain 

41 

0 

E,  Galdon,  Observatory  of  Ebro 

S.W.G 

Athene,  Greece 

38 

23 

D.  Matsoukas,  U.  of  Athens 

S.W.G 

Haifa,  Iarael 

33 

35 

J.  Mass.,  Haifa  Radio  Observatory 

S.W.G 

Addle  Ababa,  Ethiopia 

9 

38 

P.  Gouln,  U.  of  Addis  Ababa 

S.W 

Accra,  Ghana 

6 

0 

J.  Koster,  U.  of  Ghana 

S,  G 

'Nairobi,  Kenya 

-i 

37 

R.  Kelleher,  U.  of  Nairobi 

S.  G 

Dar  Es  Salaam 

-7 

39 

D.  Osborne,  U.  of  Dar  Es  Salaam 

S.W 

|  A sian  -  Pacific  Sector 

East 

Principal  Investigator 

Data 

Station  Name 

Latitude 

Longitude 

and  Affiliation 

Available 

Tehran,  Iran 

36 

51 

K.iAfshar,  U.  of  Tehran 

S 

Tokyo,  Japan 

35 

139 

Y.  Nakata,  Radio  Research  Labs 

G 

New  Dehli,  India 

28 

77 

A.  P.  Mitra,  National  Physical  Lab 

S.W 

Taipei,  Taiwan 

25 

121 

K.  Pal,  National  Taiwan  U.  ' 

S,  G 

Ahmedabad,  India 

23 

72  > 

G.  Rastogi,  Physical  Research  Lab 

S 

.Calcutta,  India 

23 

89 

S.  Basu,  U,  of  Calcutta 

S 

Hong  Kong 

22 

114 

G.  Walker,  U.  of  Hong  Kong 

S,  G 

Honolulu,  Hawaii 

21 

202 

P.  Yuen,  U.  of  Hawaii 

G 

Hyderabad,  India 

17 

78 

E.B.  Rao,  Defense  Electronics 

S 

Research  Lab 

Bangkok,  Thailand 

13 

100 

C.  Rufenach,  National  Oceanographic 

S.W 

and  Atmospheric  Administration 

Singapore 

1 

104 

E.  Golton,  Radio  Research  Station 

S 

Slough,  England 

Brisbane*  Queensland 

-27 

153 

G.  Bowman,  V.  of  Queensland 

S.W 

Armidale,  N.S.  W.  Australia 

-30 

152 

F;  Hlbberd,  U,  of  New  England 

S,  G  I 

Adelaide,  South  Australia 

-35 

138 

B.  Briggs,  U,  of  Adelaide 

? 

Auckland,  New  Zea'jnd* 

-37 

175 

J,  Titheridge,  U.  of  .Auckland 

S,  G 

Sydney,  Australia  . 

-37 

110 

G.  Munro,  U.  of  Sydney 

s 

Melbourne,  Australia 

-38 

145 

E.  Essex,  LaTrobe  University 

s,  G 

25 


Table  4-2.  (Continued) 


East 

9tation  Name  Latitude  Longitude 

Principal  Investigator  Data 

and  AfflUatlon  Available 

WeUlngton,  New  Zealand  -41  174 

Christchurch,  New  Zealand  -43  173 

Invercargill,  New  Zealand  -46  1 68 

Campbell  bland,  New  Zealand  -S3  169 

J.  Tltheridge,  U.  of  Auckland  S,  G 

G.  Stuart,  Dept,  of  Scientific  S,  7 

and  Industrial  Research 

J.  Tltheridge,  U.  of  Auckland  S,  G 

J.  Mawdsley,  Dominion  Physical  S, 

Labcratory,  DStR 

Key: 

S  indicate*  data  available  from  low-crbltlng  aatellitea 
G  indicates  data  available  from  geostationary  satellites 
W  Indicates  data  or  information  may  be  available  from  the  World  Data  Center 
*  indicates  the  principal  investigator  for  this  station  may  hove  limited  data 
from  other  stations  not  listed. 
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5.  A  Numerical  Model  of  TEC  Over  Europe 
for  Sunspot  Minimum  Conditions 

JA  Klobuchar 


5-1  INTRODUCTION 

Results  of  two  studies  of  the  TEC  of  the  ionosphere  over  Europe  have  been 
given  by  Amayenc  (1971).  They  presented  contours  of  equal  TEC  values  as  a 
function  of  geographic  latitude  and  local  time  over  the  40  to  72  deg  geographic 
latitude  range  over  Europe  for  nine  seasons  in  the  1964  to  1967  period.  The  TEC 
data  shown  in  their  iso-contour  maps  was  taken  from  two  stations  in  Europe, 
from  Kiruna,  Sweden,  and  from  Val  Joyeux,  near  Paris,  France.  Original 
glossy  prints  of  the  iso-contours  in.the  Amayenc  et  al.  paper  were  kindly  provided 
by  F.  Bertin  for  use  in  this  report.  TEC  hourly  values  were  scaled  from  these 
glossy  prints  at  5  deg  geographic  latitude  intervals  and  a  Fourier  time  series 
expansion  of  harmonic  number  4  was  made  to  each  latitude  set.  A  least  squares 
third  degree  polynomial  then  was  fitted  separately  to  each  Fourier  coefficient  over 
the  latitude  range  from  40  to  70  deg.  The  resultant  coefficients  equal  36  in  number, 
one  set  of  4  polynomial  coefficients  to  represent  the  latitude  dependence  of  each  of 
the  9  Fourier  terms.  Thus,  a  set  of  36  numbers  specify  the  seasonal  mean  TEC 
behavior  over  the  40  to  70  deg  geographic  latitude  range.  All  that  is  necessary  to 
obtain  a  value  of  TEC  is  to  use  the  coefficients  for  the  season  desired,  and  to 
specify  a  latitude  and  local  time.  The;model  is  in  the  following  form: 

4 

TEC  =  DC  +  2  £  C£  cos  (it  -  211*^/24) 
i=l 

where 

DC,  C  and  <f>  =  KQ  +  J  Kj  (lat)^  for  the  appropriate  DC,  Cj^  or  ^  term. 

3=1 


fteceiMW8  blank 
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The  coefficient  sets  for  the  seasons  are  given  in  Table  5-1.  Table  5-2  contains 
a  list  of  the  months  in  each  season  along  with  the  mean  observed  10.  7  cm  solar 
flux  over  each  season.  Note  that  themean  10.7  cm  solar  flux  was  nearly  the  same 
for  all  4  seasons. 

In  order  to  determine  how  well  the  numerical  coefficients  represent  the  actual 
data,  the  TEC  values  obtained  from  the  model  were  subtracted  from  the  original 
contour  values  at  5  deg  latitude  intervals.  The  resultant  errors  are  shown  for 
the  four  seasons  in  Figure  5-1.  Differences  less  than  1  x  1016  were  left  blank  as 
this  is  probably  the  experimental  accuracy  of  the  original  data.  Note  that  the 
model  fit  is  good  indeed,  as  evidenced  by  the  large  blank  areas  in  the  figure. 

Figure  5-2  is  an  iso-contour  of  TEC  as  constructed  by  the  numerical  model  for 
season  1  of  the  four  seasons.  It  compares  well  with  the  original  iso-contour  for 
this  season  which  is  shown  in  Figure  5-2. 

This  model  shouldr.have  use  in  determining  average  TEC  values  for  a  given 
season  for  time  delays  for  satellite  navigation  and  satellite  detection  radars 
operating  in  the  VHF  to  UHF  bands.  It  is  not  intended  for  use  in  a  day  to  day 
operational  mode,  but  as  an  average  background  model.  It  will  be  particularly 
useful  as  a  representative  TEC  background  model  for  the  approaching  1974  to  1975 
solar  minimum  conditions  for  system  design  studies  for  the  European  sector.  Its 
use  in  polar  or  equatorial  latitudes  or  at  other  than  European  longitudes  is  not 
recommended. 


Table  5-1.  Seasonal  Coefficients  Sets 


season  »ou.  tonc  se 

t  CONSTANT  -.3662*02 
t  *  »Z562*01 
1  X**2  -.6902-01 

1  x*f  3  .3062-03 

2  CONSTANT  -.2392*02 
2  X  .1972*01 
2  X**2  -.3662-01 

2  X**3  .2252-93 

3  CONSTANT  .2092*02 
3  X  -.6622*00 
3  x**2  .HIE-At 
3  X**3  -.SIX-1* 
A  CONSTANT  -.3062*02 
A  X  ,2292*01 
A  X**2  -.6512-01 
A  1X”3  .2012-03 


Cl  C2  C3  CA 

.3052*02  -.0712*01  .5912*01  .6672*01 

-.1622*01  .5322*00  •< 365E+0O  -.AttE*00 

•2992-61  *012-02  .2692-02  .9032-02 

-.I0AE-03  .550E-0A  -.A23E-0A  -.5502-06 

.2002*02  .1032*02  .3602*01  .6232*01 

-.1052*01  -.0522*00  -.1902*00  -.6132*00 
.1922-01  .1302*01)  .3602-02  .0002-02 

-.1102-03  -.5002-06  -.2132-06  -.5102-06 
.1602*02  -.2012*02  -.0352*01  -.2522*01 
-.0602*00  .1072*01  .6592*00  .2092*00 

.9502-02  -.3102-01  -.9252-02  -.6152-02 
-. 6592-06  .1002-03  .0072-06  .2532-06 

.1502*02  -.1072-01  3532*01  -.6192*01 

-.7502*00  -.1722-01  .1052*00  .2352*00 

.1362-01  .2052-03  -.3062-02  -.6362-02 

-.0102-66  -.6652-05  .1622-06  .2032-06 


FHU  FHI2  FHI3  FHI6 

.3372*02  .7172*02  -.1392*06  -.1522*02 

-.1<3E*01  -.6632*01  .0012*02  .502E*20 

.2672-01  .9022-01  -.150E«01  .1002-01 

-.1592-03  -.5022-03  .9102-02  -.3032-03 

.6572*02  -.9012*01  .2292*03  .1062*03 

-.1032*01  .7252*00  -.1602*02  -.003E*01 

.3032-01  -.3622-02  .2062*00  .1352*00 

-.2612-03  -.2222-06  -.190E-02  -.7532-03 
-.7592*02'  .5772*02  -.6292*03  .6022*03 

.5562*01  -.2032*01  .2C62*02  -.202E*02 

-.1032*00  .590E-01  -.5162*00  .6602*00 

.0022-03  -.6052-03  .3262-02  -.2702-02 

-.6572*02  -.5702*03  -.7102*03  .6012*03 

.3(02*01  .3202*02  .6102*02  -.2262*02 

-.7302-01  -.5552*00  -.7(02*00  .6102*00 

.6662-03  .3132-02  .6722-02  -.2522-02 


29 


Table  5-2.  Seasonal  Solar  Flux 


Season,. 

Month,, 

2800  MHz 
Observed  Flux 

Mean  Seasonal 
Observed  Flux 

November,  1964  . 

72.8 

I 

December,  1964 

77.5 

75.9 

January,  1965 

77.5 

February,  1965 

74.6 

II 

March,  1965 

April,  1965 

73.8 

73.4 

May,  1965 

77.9 

III 

June,  1965 

77.0 

76.4 

July,  1965 

74.3 

August,  1965 

74.8 

- 

IV 

September,  1965 
October,  1965 

76.3 

79.6 

76.9 

Winter.  1964-1965 


Latitude  0  1 

2 

3 

4 

5 

6  7 

8 

9  10 

70 

-1  -1 

1 

65 

1 

1 

1 

60 

1 

1 

-1  -1 

1 

55  -1 

50  -1 

45 

40 

1 

1 

-1 

-1 

1 

Spring.  1965 

Hour 

Latitude  0  1 

2 

3 

4 

5 

6  7 

8 

9  10 

70  1 

65 

60 

55 

50 

45 

40 

Summer.  1965 

Hour 

Latitude  0  1 

2 

3 

4 

5 

6  7 

8 

9  10 

70 

65 

1 

-1 

60 

55 

-1 

-1 

-1' 

50  -1 

45  -1 

1 

1 

1 

-1 

-1 

40  -1 

1 

-1 

1 

Autumn,  1965 
Latitude  0  1 

2 

3 

4 

5 

6  7 

8 

9  10 

70 

65 

60 

1 

1 

1 

55  -1  -1 

50 

45 

40 

1 

1 

1 

-1 

11  12  13  14  15  1<  IT  18  18  20  21  22  23 

-1  -1  1.  -1  -1  -1 
-1  1 

11  11  -1-1 

.1  -1  1 

-1  1  1 


11  12  13  14  15  16  17  18  19  20  21  22  23 

-1 

-1  -1  -1  -1 
1 

1  1 

-1 


11  12  13  14  15  16  17  18  19  20  21  22  23 


-1 


1111 
-2  -1 


11  12  13  14  15  16  17  18  18  20  21  22  23 


-1-11  1 
-1-1  111111  -1 
!  -1  -1  -1  -1 

1 


Figure  5-1.  Difference  Between  Model  and  Actual  TEC  (Units  of  10*®  el/m^) 
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WINTER  NM  •  IMS 


LOCAL  SOLAN  TIME 


Figure  5-2a.  Iso-Contour  of  TEC 
Northern  Europe  as  Reconstructed  by 
the  Numerical  Model,  Winter 
1964-65 


WINTER  IM4-IM5 


Figure  5-2b.  Original  Iso-Contour  of 
TEC  Over  Northern  Europe,  Winter 
1964-65 
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6.  A  Comparison  of  TEC  Obtained  From 
Thomson  Scatter  and  From  Faraday  Rotation 

J.A.  Klobachar  and  S.  Basu 


6-1  A  COMPARISON  OF  TEC  COMBINED  FROM  THOMSON 

SCATTER  AND  FROM  FARADAY  ROTATION 

Continuous  measurements  of  TEC  have.been  made  for  several  years  at 
Sagamore  Hill,  Hamilton,  Massachusetts  by  observing  the  Faraday  polarization 
twist  of  VHF  radio  waves  transmitted  from  the  geostationary  satellite  ATS- 3.  The 
ATS-3  satellite  has  been  moved  in  longitude  from  time  to  time,  but  has  spent  much 
of  its  time  nearly  south  of  Hamilton,  Mass,  at  an  elevation  of  approximately  40  deg. 
The  latitude  below  the  point  where  the  ray  from  Hamilton,  Mass,  to  the  satellite 
crosses  the  region  of  maximum  electron  density’is  approximately  39  deg  North 
geographic. 

Routine  measurements  of  vertical  ionospheric  electron  density  profiles  were 
made  at  Millstone  Hill,  Westford,  Mass.,  42.6  deg  North,  by  J.  V.  Evans,  Lincoln 
Laboratory,  during  a  few  days  per  month.  In  this  section  of  the  report  comparisons 
are  made  between  four  days  of  vertical  TEC  taken  by  the  Thomson  scatter  technique, 
kindly  supplied  by  J.  V.  Evans,  and  "equivalent  vertical  TEC"  obtained  by  looking 
at  the  Faraday  polarization  twist  along  the  slant  path  to  the  ATS-3  satellite.  The 
four  days  used  in  this  comparison  were  magnetically  quiet  days,  probably  typical  of 
the  four  different  seasons  during  a  sunspot  maximum  year. 

The  Thomson  scatter  data  were  usually  available  from  a  height  below  200  km  to 
1200  km,  though  during  some  nighttime  periods  the  signal -to-noise  ratio  was  too 
poor  to  obtain  data  to  1200  km  due  to  the  low  densities  at  very  great  heights.  The 
TEC  was  obtained  from  the  Thomson  scatter  data  simply  by  integrating  the  profiles 
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throughout  the  height  range  for  which  data  was  provided.  No  extrapolation  above 
1200  km  was  made,  nor  was  any  attempt  made  to  fill  in  data  below  the  minimum 
height  given  by  the  Thomson  scatter  data.  The  profiles  were  integrated  at  approxi¬ 
mate  hourly  intervals,  during  each  of  the  four  days. 

The  TEC  derived  from  the  Faraday  effect,  was  obtained  by-converting  the 
observed  Faraday  polarization  twist  data  to  equivalent  vertical  TEC  by  using  the 
longitudinal  magnetic  field  and  zenith  angle  at  a  fixed  mean  ionospheric  height  of 
350  km.  Results  of  model  calculations  have  shown  that  the  Faraday, data  obtain  TEC 
out  to  heights  of  approximately  2000  km.  Thus,  the  Faraday  TEC  should  be  greater 
than  the  Thomson  scatter  TEC,  because  it  is  integrated  to  approximately  2000  km, 
whereas  the  Thomson  scatter  TEC  is  integrated -to- a -maximum  of  1200  km.  Also, 
there  is  an  approximate  3  deg  latitude  difference  between  the  locations  of  the  two 
measurements,  with  the  Faraday  TEC  being  the  farthest  south.  In  this  northern 
mid-latitude  region  gradients  of  electron  density  are  usually  lower  at  higher  latitudes. 
These  two  factors  should  combine  to  make  the  Faraday  TEC  higher  than  the  TEC 
obtained:  from  the  Thomson  scater. 

Results  of  the  four  days  of  comparisons  are  shown  in  Figure  6-1.  The  agree¬ 
ment  between  the  two  sets  of  data  on  all  four  days  is  excellent.  The  largest  differ¬ 
ence  between  the  two  data  sets  occurs  during  the  daytime  afternoon  periods  when  the 
latitude  gradient  of  electron  density  between  42. 6  deg  North  and  39  deg  North  is 
expected  to  be  the  largest.  No  attempt  has  been  made  to  correct  either  data  set  for 
the  effects  of  the  expected  latitude  gradient.  While  the  difference  between  the  two 
data  sets  in  the  daytime  period  is  likely  due  to  the  latitude  gradient  between  the  two 
observation  latitudes,  the  observed  difference  at  nighttime  when  the  latitude  gradient 
is  small  is  probably  due  to  the  additional  contribution  to  TEC  above  1200  km  in  the 
Faraday  data.  This  additional  contribution  also  accounts  for  part  of  the  daytime 
additional  content  observed  in  the  Faraday  data. 

The  one  nighttime  case  when  the  Thomson  scatter  TEC  was  higher  than  the 
Faraday  TEC  occurred  in  the  pre-dawn  hours  of  January  17,  1969.  The  winter 
nighttime  F  region  sometimes  has  relatively  large  gradients  having  small  geographic 
extent.  The  pre-dawn  dip  observed  in  Faraday  TEC,  while  a  small  nighttime 
increase  was  observed  in  Thomson  scatter  TEC  was  likely  the  result  of  such  a  winter 
nighttime  gradient.  Despite  the  small  differences  due  to  the  reasons  cited  above,  the 
TEC  taken  by  the  Faraday  polarization  twist  method  agrees  well  with  the  TEC 
obtained  by  integrating  electron  density  profiles  obtained  by  Thomson  scatter. 
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7-1  INTRODUCTION 

to  this  detasection  of  the  report,  typical  behavior  ot  TEC.  and.  where  available, 
slab  thickness  values  are  presented,  along  with  a  tew  comments  on  the  general 
state  of  knowledge  of  TEC  in  each  region. 


7-2  EQUATORIAL  TEC  AND  SLAB  THICKNESS 

The  TEC  oft  the  equatorial  region  generally  follows  the  same  anomalous  behavior 
as  the  density  of  the  peak  of  the  F  region,  Nmax,  that  is,  the  TEC  and=Nmax  are  not 
highest  over  the  sub-solar  latitude  (the  equator,  for  example,  at  the  equinoxes). 

Both  these  parameters  generally  are  highest  at  latitudes  displaced  on  both  sides  of 

the  dip  equator.  ,  .  , 

One  of  the  best  summaries  of  TEC  in  the  equatorial  region  was  made  by  de 

Mendonca,  et  al  (1969)  by  observing  the  Faraday  effect  from  VHF  signals  transmitted 
from  S-66  (BE-B  and  BE-C)  satellites  over  Brazil  in  1966.  Their  data  was  reduced 
in  detail  and  plotted  in  the  form  of  contour  maps  of  constant  TEC  values  over  the 
dip  latitude  0  deg  to  35  deg  South  and  0-24  hours  local  mean  time.  Figure  7-1  shows 
de  Mendonca's  latitude-time  contours  of  TEC  for  four  seasons. 

At  Hawaii,  which  is  located  on  the  northern  edge  of  the  equatorial  anomaly,  a 
several  year  history  of  continuous  TEC  data  has  been  compiled  by  monitoring 
Faraday  polarization  from  VHF  signals  transmitted  from  one  or  more  geostationary 
satellites.  Yuen  and  Roelofs  (1967)  show  much  of  this  data. 
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Figure  7-1.  Contour  Maps  of  Total  Electron  Content  for  the  Four  Seasons 

of  1966  (BE-B  and  B-C) 

Results  of  one  year's  TEC  observations  taken  at  the  Panama  Canal  Zone  are 
shown  in  Figure  7-2.  Despite  the  loss  of  much  of  the  data  due  to  numerous  equip¬ 
ment  and  personnel  problems,  during  some  months  there  was  sufficient  data  taken 
to  determine  a  useful  median  TEC  value.  Since  there  was  no  nearby  ground  based 
ionoscnde  for  finding  Nmax,  monthly  median  values  of  slab  thickness  were  computed 
using  the  median  TEC  values  and  predicted  values  of  Nmax  from  the.ESSA  numerical 
coefficient  model.  These  slab  thickness  values  are  shown  in  Figure  7-3.  Since 
these  curves  were  constructed  using  monthly  median  values  of  TEC  obtained  from 
only  a  few  days  of  data  per  month  along  with  numerically  predicted  values  of  N  .  . 
they  are  of  low  accuracy. 

Other  studies  of  TEC  in  the  equatorial  region,  such  as  those  by  Basu  and  das 
Gupta  (1967),  Rufenach,  et  al  (1968),  Hunterv  pt^al  (1965),  and  Blumle  (1962)  were 
of  short  time  duration  or  were  of  more  specialized  interest  in  a  particular  feature 
of  the  equatorial  TEC.  Where  available  from  published  papers,  values  of  equatorial 
slab  thickness  have  been  included  in  a  previous  section  of  this  report. 

Figure  7-4  shows  the  TEC  behavior  as  observed  from  Kingston,  Jamaica,  West 
Indies  for  March  1972.  Faraday  rotation  observations  from  the  VHF  signals  from 
ATS-3  are  being  made  routinely  from  this  location,  and  the  data  shown  in  Figure  7-4, 
kindly  provided  by  P.  Chin,  University  of  the  West  Indies,  represents  the  first  data 
in  final  TEC  form  from  this  station.  Additional  continuous,  reduced  TEC  data  from 
this  stajidn  should  soon  become  available. 

The  TEC  contours  taken  from  Brazil  and  shown  in  Figure  7-1  probably  best 
summarize  TEC  behavior  in  the  equatorial  regions.  Near  specific  low  latitude 
stations,  such  as  Hawaii,  more  complete  TEC  data  is  available.  More  continuous 
TEC  data  should  soon  be  available  from  Jamaica,  Peru,  Puerto  Rico,  Brazil, 
and  perhaps  other  low  sites  to  eventually  allow  a  more  complete  picture  of  equatorial 
TEC  behavior  to  be  made. 


LOCAL  MEAN  TIME 


Figure  7-2.  Total  Electron  Content  for  1969  From  the  Panama  Canal  Zone 

7-3  MID-LATITUDE  TEC  AND  SLAB  THICKNESS  DATA 

The  TEC  of  the  mid-latitude  ionosphere  is  becomming  fairly  well  known,  at 
least  over  North  America,  where  several  stations  have  been  making  continuous  TEC 
measurements  for  at  least  5  years.  A  first  order  model  of  mid-latitude  slab  thick¬ 
ness  has  been  discussed  in  this  report  as  have  several  TEC  models  or  complete 
vertical  profile  models  from  which  TEC  may  be  found.  While  these  models  may  be 
available  to  system  design  engineers  for  correction  for  the  effects  of  TEC  on  their 
system,  these  models  say  little  or  nothing  about  the  statistical  behavior  of  TEC, 
that  is,  its  day  to  day  variability.  For  this  reason  the  following  tabulated  values  of 
TEC  are  presented  in  this  data  section.  Tables  7-1  through  7-6  givedhe  mean, 
standard  deviation,  minimum  and  maximum  values,  median,  upper  and  lower 
quartiles  and  upper  and  lower  deciles  of  each  hour  of  TEC  and  of  equivalent  slab 
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Table  7-1.  Hourly  Statistical  Values  of  Total  Electron  Content  and  Equivalent  Slab 
Thickness  From  Hamilton,  Massachusetts,  for  January-February  1968 
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Table  7-2.  Hourly  Statistical  Values  of  Total  Electron  Content  and  Equivalent  Slab 
Thickness  From  HamUton,  Massachusetts,  for  March-April  1968 
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Table  7-3.  Hourly  Statistical  Values  of  Total  Electron  Content  and  Equivalent  Slab 
Thickness  From  Hamilton,  Massachusetts,  for  May-June  1968 
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Table  7-4.  Hourly  Statistical  Values  of  Total  Electron  Content  and  Equivalent  Slab 
Thickness  From  Hamilton,  Massachusetts,  for  July-August  1968 
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Table  7  -6.  Hourly  Statistical  Values  of  Total  Electron  Content  and  Equivalent  Slab 
Thickness  From  Hamilton,  Massachusetts,  for  November-December  1968 
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Figure  7-3.  Monthly  Median  Equivalent  Slab  Thickness  for  1969  From  the 
Panama  Canal  Zone 
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Figure  7-4.  Total  Electron  Content 
From  Kingston,  Jamaica,  West 
Indies,  March  1972 


thickness  separately  for  each  month  of 
the  year  1968  for  data  taken  at  Hamilton, 
Massachusetts.  With  these  tabulated 
values,  which  are  typical  of  mid-latitudes, 
a  design  engineer  can  determine,  in  addi¬ 
tion  to  the  average  behavior,  the  complete 
statistical  behavior  of  the  mid-latitude 
TEC. 

In  order  to  show  the  relative  good 
agreement  between,  slab  thickness  values 
taken  at  widely  separated  mid-latitude 
stations,  a  series  of  graphs  of  slab  thick¬ 
ness  at  pairs  of' stations  is  shown  in 
Figures  7-5  through  tf-8.  All  these 
stations  are  in  the  northern  mid-latitude 
region.  These  curves  show  that  a  TEC 
model  based  upon  slab  thickness  can 
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01  12  -  II 

LOCAL  MEAN  TIME 
MEDIAN  SLAI  THICKNESS 


Figure  7-5.  A  Comparison  of  Seasonal 
Median  Values  of  Slab  Thickness  at 
Florence,  Italy,  and  Hamilton, 
Massachusetts 


Figure  7-6.  A  Comparison  of  Seasonal 
Median  Values  of  Slab  Thickness  at 
Florence,  Italy,  and  Hamilton, 
Massachusetts 


Figure  7 -7a.  A  Comparison  of  Monthly 
Median  Slab  Thickness  at  Aberystwyth, 
Wales,  and  Hamilton,  Massachusetts 


LOCAL  MEAN  TIME 
MEDIAN  SLAB  THICKNESS 


Figure  7-7b.  A  Comparison  of  Monthly 
Median  Slab  Thickness  at  Stanford, 
California,  and  Hamilton,  Massachusetts 

apply  over  this  general  region.  The 
Stanford  data  was  supplied  through  the 
courtesy  of  A.  V.  daRosa,  the  data  from 
Aberystwyth,  Wales  by  L,  Kersley  and 
the  data  from  Florence,  Italy  is  from 
P.  F.  Checcacci  and  M.  T.  deGiorgio  (all 
private  communications). 

To  give  a  picture  of  the  average 
latitude  and  temporal  behavior  of  the  mid¬ 
latitude  TEC  for  various  seasons  Figures 
7-9  and  7-10  show  TEXD  values  taken 
from  Hamilton,  Massachusetts  from  the 
VHF  signals  from  the  S-66  (BE-B)  satellite. 
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Figure  7-8.  A  Comparison  of  Monthly  Median  Slab  Thickness  at  Stanford, 
California,  and  Hamilton,  Massachusetts 

Using  the  same  satellite,  observations  of  TEC  have  been  made  from  Athens,  Greece 
and  are  shown  in  Figure  7-11.  These  are  kindly  provided  by  D,  Matsoukas  (private 
communication).  Figures  7-12  and  7-13  show  the  behavior  of  TEC  over  northern 
Europe  for  eight  seasons.  They  are  the  result  of  combining  data  taken  by  Liszka 
(1967)  with  those  taken  at  Val  Joyeux,  near  Paris,  France,  Amayenc  (1971).  These 
contours  were  kindly  supplied  by  F.  Bertin  for  inclusion  in  this  report.  The  first 
four  of  these  contours  were  used  to  make  the  numerical  model  of  TEC  over  Europe 
for  sunspot  minimum  conditions,  the  details  of  which  are  in  another  section  of  this 
report. 

While  much  work  on  TEC'behavior  has  been  done  at  a  few  other  mid-latitude 
locations,  such  as  in  New  Zealand,  by  Titheridge,  and  in  Japan  by  Nakata,  the 
stress  on  data  presented  here  is  on  the  Atlantic  and  European  sectors  because  of 
NATO's  particular  interests  in  this  region.  The  statistics  given  in  Tables  7-1 
through  7-6,  along  with  the  many  TEC  contours  shown  in  Figures  7-9  through  7-13, 
depict  fairly  well  the  behavior  of  TEC  in  this  region. 
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Figure  7-9i  Iso-Contours  of  Seasonal  Median  Total  Electron  Content  Taken  at 
Sagamore  Hill,  Hamilton,  Massachusetts 


7-4.  THE  HIGH  LATITUDE  REGION 

TEC  measurements  at  the  high  latitudes  have  been  made  by  only  a  few  workers. 
Some  of  the  important,  studies  and  their  data  are  summarized  here.  TEC  measure¬ 
ments  using  the  S-66  (BE-B),  1000  km  height,  satellite  were  made  from  Kiruna, 
Sweden  (67  deg  North,  20  deg  East)  geographic  coordinates  by  Liszka  (1967),  A 
trough  in  the  latitude  distribution  of  TEC  was  sometimes  observed  on  individual 
passes,  particularly  during  the  winter  nighttime  when  it  was  always  observed. 
Because  the  satellite  is  generally  viewed  at  an  oblique  angle,  sharp  gradients  in 
TEC  are  smoothed  by  this  measuring  technique.  Liszka  constructed  iso-contour 
maps  of  TEC  versus  latitude  and  local  time  for  all  hours  for  9  three-month  seasons 
over  the  latitude  range  from  55  to  72  deg  North  over  the  European  sector.  These 
contours  are  shown  in  the  mid-latitude  portion  of  the  data  section  as  they  were 
attached  to  the  mid -latitude  TEC  contours  taken  near  Paris  by  Amayenc,  et  al.  The 
trough  in  TEC  is  clearly  seen  in  both  the  winter  nighttime  contours  of  the  winter 
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Figure  7-10.  Iso-Contours  of  Seasonal  Median  Total  Electron  Content  Taken  at 
Sagamore  Hill,  Hamilton,  Massachusetts 

1964  to  1965  and  the  winter  of  1965  to  1966.  Though  there  is  some  evidence  of  the 
trough  during  the  equinox  periods,  the  rapid  seasonal  changes  which,  occur  at  these 
high  geogr^iiqal  latitudes  were  responsible  for  averaging  both  day  and  night  values 
into  th<?  equinox  contours,  thus  smoothing  out  any  nighttime  trough.  Contours  of  the 
equivalent  slab  thickness  parameter,  over  a  10  deg  latitude  range,  were  constructed 
by  Liszka  for  four  seasons  by  combining  his  TEC  observations  with  foF2  observa¬ 
tions  from  six  Scandinavian  ionosondes.  These  results  are  shown  in- Figure  7-14. 

In  the  winter  two  distinct  maxima  in  slab  thickness  are  clearly  visible,  one  at 
21  hours  and  the  second  at  03  hours  local  time.  Both  occur  at  65  deg  Northgeo- 
graphic  latitude.  These  maxima  occur  when  the  trough  is  above  the  Kiruna  station. 

Bratteng  and  Frihagen  (1969)  reported  the  results  of  TEC  measurements  made 
at  Spitzbergen  (78  deg  North,  14  deg  East).  They  constructed  iso-contours  of  TEC 
versus  local  time  for  the  Spring  and  Autumn  1965  periods.  Their  contours  are 
reproduced  in  Figure  7-15.  They  noted  the  large  0900  LMT  maximum  in  TEC  at 
78  deg  North,  but  little  evidence  of'a>trough  during  the  seasons  shown  is  noticeable. 
Mikkelsen  (1971)  combined  observations  of  TEC  taken  at  St.  John’s,  Newfoundland 
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Figure  7-11.  Iso-Contours  of  Seasonal  Median  Total  Electron  Content  Taken 
at  Athens,  Greece 

and  Narssarssuaq,  Greenland  during  July  1969  using  the,S-66  beacon  satellite 
transmissions.  He  found  a  knee  in  TEC  which  moved  southward  with  increasing 
magnetic  activity,  but  not  very  good  agreement  between  the  position  of  the  knee  and 
that  of  the  southern  edge  of  the  scintillation  boundary.  No  mean  data  was  shown  by 
Mikkelsen  either  for  TEC  versus  latitude  or  time. 

Rai  and  Hook  (1967)  made  TEC  observations  near  College,  Alaska  during  the 
winter  of  1963-1964.  They  found  a  large  secondary  nighttime  maximum  in  TEC 
near  local  midnight  between  65  deg  and  75  deg  North  geographic  latitude.  At  60  deg 
North  there  was  only  a  slight  nighttime  increase.  During  magnetically  disturbed 
conditions,  they  found  large  positive  gradients  towards  the  North.  They  reported 
no  measurements  of  slab  thickness. 

Currently,  continuous  measurements  of  TEC  are  being  taken  at  several  high- 
latitude  locations  using  the  VHF  signals  transmitted  from  the  geostationary  satellite 
ATS-3.  Stations  where  reduced  data  will  soon  be  available  include  Narssarssuaq, 
Greenland,  Goose  Bay,  Labrador,  and  College,  Alaska.  Some  three  months  of 
data  from  Thule,  Greenland  have  been  reduced  and  are  outlined  in  the  following 
section. 


Figure  7-14.  Contours  of  Equivalent  JSlab 
Thickness  Over  Northern  Europe  -  1965 
(after  Liszka,  1967) 
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Figure  7-15a.  Contour  Map  of  Electron  Content  in  a  Geographic 
Latitude-LMT  Plane  for  the  Period  15/3-15/4/66  (After  Bratteng 
and  Frihagen,  1969) 
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Figure  7-15b,  Contour  Map  of  Electron  Content  in  a  Geo¬ 
graphic  Latitude-LMT  Plane  for  the  Period  15/10-15/11/65 
(After  Bratteng  and  Frihagen,  1969) 
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8.  Low  Elevation  Angle  Measurements  of  Total 

Electron  Content  Taken  From 
Thule,  Greenland 

M.^MtfldHIo  and  J.A.  Klobuchar 


8-1  INTRODUCTION 

Continuous  measurements  of  the  Faraday  polarization  twist  of  VHF  radio  waves 
from  the  geostationary  satellite  ATS-3  have  been  made  from  April  1971,  at  Thule, 
Greenland  (75. 5  deg  North,  68. 7  deg  West).  This  set  of  data  represents  the  first 
attempt  made  from  such  a  high  latitude  site  to  continuously  monitor  the  diurnal 
behavior  of  the  ionospheric  total  electron  content  (TEC).  During  this  period  the 
ATS-3  satellite  was  close  to  the  Thule  meridian  (70  deg  West)  at  a  ground  elevation 
angle  of  approximately  5  deg.  Figure  0-1  shows  the  sub-ionospheric  geographic 
latitude  versus  ionospheric  height  together  with  the  L  value  at  the  Faraday  rotation 
factor  M  (BcostfsecX)  along  the  ray  path  to  ATS-3.  Note  that  the  ionospheric  station 
at  Narssarssuaq,  Greenland  (61. 2  deg  North  45. 4  deg  West),  is  located  at  the  same 
geographic  latitude  as  the  400  km  sub -ionospheric  coordinates.  The  L  shell  at  this 
point  is  11.6  which  corresponds  to  an  invariant  latitude  (a)  of  73  deg. 

The  values  of  polarization  twist  (Q)  recorded  at'  Thule  have  bden  converted  to 

"equivalent  vertical"  electron  contents  (N^,)  using  Q  =  —  MN,p  with  the  mean  M 

factor  computed  for  an  ionospheric  height  of  400  km.  It  should  be  noted  that  over 
the  height  range  300  to  500  km  M  changes  by  only  7  percent,  and  thus  the  decision 
at  which  height  M  should  be  chosen  has  more  of  an  effect  on  determining  the  sub- 
ionospheric  coordinates  that  it  does  on  the  derived  values  of  Nq-.  The  polarimeter 
in  use  at  Thule  has  been  calibrated  orily  tentatively  for  absolute  polarization  reference, 
and  therefore  the  absolute  scale  for  the  TEC  curves  was  arbitrarily  chosen  to  give 
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reasonable  values  of  slab-thickness 
(  t  =  N<p  /  Nmax).  It  should  be  empha¬ 
sized,  however,  that  the  data  are  con¬ 
tinuous  over  an  approximate  three  month 
period  and  thus  day-to-day  and  seasonal 
differences  may  confidently  be  compared 
on  a  relative  basis. 

Before  presenting  our  preliminary 
results,  we  would  like  to  comment  briefly 
on  the  appropriateness  of  constructing 
"equivalent  vertical"  TEC  values  using 
the  Faraday  rotation  technique  with  a 
satellite  only  5  deg  off  the  horizon. 

Model  calculations  were  performed  for 
day  and  nighttime  equinox  conditions 
using  the  geometry  depicted  in  Figure  8-1. 
The  electron  density  versus  height  and 
latitude  were  taken  from  Chan  and  Colin  (1969)  for  the  topside  with  extrapolated 
values  for  the  bottomside.  When  the  amounts  of  rotation  computed  along  the  slant 
path  through  the  ionosphere  were  converted  to  equivalent  vertical  N«p  values  using 
M  at  400  km,  the  daytime  results  agreed  to  within  3  percent  and  the  nighttime  to 
within  7  percent  of  the  actual  vertical  integrated  profiles  passing  through  the  400  km 
sub-ionospheric  point.  Several  of  the  factors  which  contribute  to  this  good  agree¬ 
ment  are: 

(1)  Since  TvT  and  Ne  decrease  both  with  height  and  to  the  south,  most  of  the 
rotation  occurs  where  the  ray  path  intersects  the  F2  peak. 

(2)  The  latitudinal  Ne  gradients  southwards  towards  the  trough  are  not  too  great 
at  ionospheric  heights  near 

(3)  While  the  elevation  angle  to  ATS-3  is  only  5  deg  at  the  ground  (sec  X  -  11 .  V, 
at  400  km  along  the  ray.  path  it  is  20  deg  (sec  X  -  3.0). 


Figure  8-1.  Sub-Ionospheric  Latitude 
and  M  Factor  vs  Vertical  Height  From 
Thule  Viewing  ATS-3  Along  the  Station 
Meridian 


8-2  DESCRIPTION  OF  THE  DATA 

In  Figure  8-2  data  from  the  first  three  months  of  operation  are  given  in  a  mass- 
plot  format.  The  tendency  for  the  TEC  to  have  preferred  values  occurs  because  the 
received  polarization  is  a  combination  of  a  direct  ray  from  the  satellite  and  a  ground 
reflected  ray  of  constant  horizontal  polarization.  The  contamination  from  the  ground 
reflected  ray  is  due  to  the  low  elevation  angle  at  which  these  observations  were  made 
using  modest  gain  yagi  antennas.  Despite  this  problem,  some  of  the  features  clearly 
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Figure  8-2.  Total  Electron  Content 
From  ATS-3,  Thule,  Greenland,  1971. 
Monthly  overplots  of,  all  data 


visible  in  Figure  8-2  include  a  diurnal 
effect  with  a  broad  afternoon  maximum, 
a  day-to-day  variability;  and  also  a  great 
deal  of  structure  at  virtually  all  local 
times.  The  monthly  median  curves  (ftp) 
are  presented  in  Figure  8-3.  One  can 
see  that  a  seasonal  trend  is  apparent  in 
that  at  all  local  times  "ftp  (June)  >  ftp. 
(May)  >  ftp.  (April)  with  the  nighttime 
values  showing  the  greatest  change.  At 
mid-latitudes,  on  the  other  hand,  the 
seasonal  variation  during  the  daytime 
hours  is  quite  different,  that  is,  daytime 
ftp.  values  in  June: are  significantly  lower 
than  the  corresponding  ftp,  values  in 
Aprii.  Figure  8-3  also  shows  that  the 
overall  diurnal  range,  that  is,  ftp(max)  - 
ftp,(min),  also  varies  with  season.  During 
April,  the  range  is  12  x  10*6  el/m^  while 


during  May  and  June  it  reduces  to  10  and 
8  x  1016  el/m2,  respectively.  Finally,  the  small  arrows  in  Figure  8-3  give  the 
times  of  ground  sunrise  and  sunset  at  61  deg  North.  It  appears  that  as  the  daylength 
increases  with  season;  the  slopes  (dNq>/dt)  at  sunrise  and  sunset  both  decrease. 


LOCAL  TIMC 


Figure  8-3.  Total  Electron  Content 
From  ATS-.3,  Thule,  Greenland,  1971. 
Monthly  medians 


8-3.  THE  INFLUENCE  OF  GEOMAGNETIC 
DISTURBANCES 

In  order  to  investigate  the  effects  of 
magnetic  activity,  we  concentrated  on  the 
April  9  to  30  period.  During  this  time, 
two  moderately  severe  sudden  commence¬ 
ment  storms  occurred  and  in  both  cases 
there  were  recognizable  departures  in 
Nt  from  monthly  median  conditions. 

From  Figures  8-4a  and  b  one  can  see 


that  both  of  the  "ionospheric  storms"  took  the  form  of  an  initial  positive  phase 
followed  by  a  longer  negative  phase.  This  type  of  storm-time  (Dst)  behavior  is 
quite  similar  to  that  normally  found  at  mid-latitudes  [Hibberd  and  Ross,  (1967);  and 
Mendillo  (1971)]  .  With  only  two  storm  periods  available,  however,  it  is  not 
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Figure  8-4.  Total  Electron  Content  from  ATS-3,  Thule*  Greenland,  During  One 
Quiet  and  Two  Disturbed  Periods  in  April  1971 

possible  to  discuss  local  time  (SD)  effects.  The  SC  for  the  first  storm  occurred  at 
00:29  LT  on  9  April  and  the  enhancement  in  N<p  occurred  during  the  broad  minimum 
usually  found  near  dawn.  The  second  SC  occurred  at  07:43  LT  on  the  14th  and  the 
positive  phase  maximum  in  NT  occurred  near  local  noon.  Thus,  in  both  cases  the 
time-lag  from  SC  to  ANT(max)  was  3  to  5  hours.  Figure  8-4c  gives  an  example 
of  some  nighttime  increases  in  NT  during  a  period  when  no  storm  SC  was  reported 
but  when  global  conditions  were  generally  disturbed.  During  the  6  hour  period 
countered  on  00  UT  of  22  April,  Kp  was  6+  and  5".  We  agree  with  the  suggestions 
made  by-Rai  and  Hook  (1967)  and  Bradbury  et  al  (1968)  that  the  precipitation  of  low 
energy  electrons  is  the  probable  cause  of  nighttime  enhancements  in  TEC  at  auroral 
latitudes. 

In  order  to  present  in  the  greatest  detail  available  a  few  examples  of  the  diurnal 
N.p  curves  measured  from  Thule,  we  have  prepared  in  Figure  3-5  tracings  of  the 

actual  polarimeter  data  taken  on  14  April  (Ap  =  39,  £  Kp  =  28+)  and  on  25  April 

(Ap  =  1,  £  Kp  =  2-).  The  storm  effects  on  the  14th  are  quite  dramatic,  while  the 

extremely  calm  (geomagnetically)  day  closely  follows  the  monthly  median  curve  in 
Figure  8-3.  It  is  interesting  to  note  that  even  on  a  quiet  day  a  certain  amount  of 
variability  is  apparent. 
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Figure  8-5.  ATS-3  From  Thule,  Greenland,  14  and  25-April  1971. 

(76. 6°N,  68. 7°W),  Sub-Ionospheric  Point  (400  km)  =  61.6°N,  75.1°W,  A  =  73°. 
Elevation  to  ATS-3  =  4. 6° 


84.  DISCUSSION 


Perhaps  the  three  most  prominent  features  evident  in  this  preliminary  set  of 
data  are 

(1)  the  solar  zenith  angle  influence  in  the  diurnal  variation, 

(2)  the  definite  response  noted  to  geomagnetic  storms,  and 

(3)  the  irregular  fluctuations  found  at  all  Hmes,  especially  the  nighttime  hours. 
There  are  at  the  present  time  no  published  reports  of  the  complete  diurnal  N-p 
behavior  at  high  latitudes  with  which  these  results  may  be  compared.  However,  for 
several  years  Liszka's  group  in  Sweden  and  Frihagen's  group  in  Norway  have  used 
low-orbit  data  to  extensively  study  total  content  at  high  latitudes.  The  Norwegian 
results  include  a  particularly  interesting  set  of  data  from  Kjeller  (60  deg  North, 

11  deg  East)  and  Spitzbergen  (78  deg  North,  14  deg  East).  NT  data  obtained  from 
the  40  and  41  MHz  beacons  on  the  satellite  S-66  as  seen  from  Kjeller  refer  to 
approximately  the  same  geographic  latitude  (60  deg)  as  the  present  Thule  data,  while 
S-66  passes  taken  near  76  deg  from  Spitzbergen  yield  Nj  data  at  approximately  the 
same  L-shell  as  the  Thule  data.  At  middle  and  low  latitudes  (that  is,  within  the 
plasmasphere)  the  geomagnetic  field  plays  such  a  dominating  role  in  determining 
the  distribution  of  ionospheric  plasma  that  one  normally  thinks  in  terms  of 
comparing  data  from  two  stations  which  have  the  same  geomagnetic  coordinates. 

The  data  taken  from  Thule,  however,  measures  electrons  situated  well  above  the 
plasmapause  (that  is,  north  of  the  main  trough)  and  thus  it  is  not  clear  whether 
geomagnetic  or  geographic  (solar)  control  has  the  greatest  influence.  Bratteng  and 
Frihagen  (1969)  presented  some  Spitzbergen  data  taken  during  the  Spring  of  1966 
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which  seems  to  indicate  that  the  overall  diurnal  range  at  76  deg  was  less  than  the 
range  reported  at  Kjeller  (60  deg)  for  the  Spring  of  1965  [Frihagen  and  Bratteng 
(1967)]  .  Considering  the  April  median  in  Figure  8-3  as  well  as  feature(l)  noted 
above,  it  appears  that  the  solar  zenith  angle  control  exerts  the  dominant  influence 
over  the  diurnal  Nj  variation-measured  from  Thule. 

In  conclusion,  we  wi3h  to  say  that  while  emphasis  has  been  placed  on  the  solar 
control  over  the  TEC  results  reported  here,  the  corpuscular  sources  of  ionization 
are  certainly  important.  This  was  briefly  discussed  in  connection  with  Figure  8-4  c 
and  it  appears  that  particle  influence  may  also  be  evident  in  Figure  8-2.  Near  1500 
LT  in  this  figure,  there  seems  to  be  an  avoidance  of  low  Nt  values.  This  time 
period  corresponds  to  approximately  2000  UT,  the  same  UT  period  when  a  maximum 
in  the  soft  electron  precipitation  is  expected  [see  Maehlum,  (1968)  and  Pike,  (1 970)] . 
Bratteng  and  Brihagen  (1969)  reported  similar  effects  during  winter  months  when 
the  possibility  of  observing  a  UT  variation  is  the  greatest.  Such  effects  will  be 
examined  in  more  detail  as  soon  as  our- Winter  data  becomes  available. 
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9.  Conclusions 

J.A.  Klobuchar 

91.  METHODS  USED,  PROBLEMS  AND  LIMITATIONS  IN  OBTAINING  TEC  DATA 

The  present  state  of  the  art  of  TEC  values  for  use  in  time  delay  corrections  for 
systems  applications  is  a  mixed  bag.  In  this  report  we  have  outlined  several  existing 
models  from  which  a  TEC  value  may  be  obtained.  We  have  tried  to  avoid  compari¬ 
sons  among  these  models  as  they  have  not  been  adequately  and  uniformly  tested 
against  real  data  in  any  comparative  way.  Our  emphasis  here  has  been  to  present 
data  from  the  three,  more  or  less,  naturally  defined  world  regions,  namely  the 
equatorial,  mid-latitude  and  polar  regions.  In  one  case  sufficient  average  data  were 
available  to  construct  a  model  of  TEC  which  is  of  use  in  the  European  sector  in  a 
.limited  way.  Also,  the  statistics  of  TEC  variability  were  given  in  the  data  section 
for  a  mid-latitude  station.  Many  maps  of  TEC  iso-contours  also  were  shown, 
particularly  for  the  European  sector..  Any  one  or  a  combination  of  the  above  models 
and/or  data  sets  should  serve  as  a  first  order  TEC  time  delay  correction  source. 

Several  problems  and  limitations  in  TEC  data  were  discussed  and  more  problems 
will  be  addressed  briefly  here.  The  basic  problems  are: 

(1)  the  lack  of  complete  measurements  of  the  many  variables  responsible  for 
the  behavior  of  the  electron  density  in  the  ionosphere;  and 

(2)  the  lack  of  sufficient  experimental  data  to  produce  an  empirical  ionospheric 
electron  density  model  of  sufficient  accuracy. 

Since  the  approach  in  (1)  follows  from  first  principles  it  is  the  more  elegant,  but 
also  far  from  complete,  solution.  Hence,  an  empirical  model  of  the  parameter  of 
direct  interest  is  the  alternate  approach,  and  all  empirical  ionospheric  models  are 
data-base  limited. 
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The  specific  TEX]  data  base  limitations,  aside  from  the  obvious  one  of  lack  of 
data  from  a  particular  geographic  region,  are  several.  TEC  from  low  orbit  satellites 
is  inadequate  because  it  can  only  determine  TEC  up  to  the  satellite  height  of  approxi¬ 
mately  1000  km,  and  then  only  when  the  satellite  is  sufficiently  above  the  station's 
horizon,  which  may  be  only  for  four  to  six  10-minute  periods  per  day.  During 
times  of  heavy  amplitude  scintillation,  the  Faraday  fading  is  obscured  and  often 
unreadable.  This  happens  particularly  during  magnetic  storms  when  the  TEC  is 
probably  most  interesting  to  study,  as  during  disturbed  periods  the  TEC  departure 
from  median  conditions  are  usually  largest.  Because  of  the  motion  of  the  satellite, 
changes  in  its  beacon  signal's  transmitted  polarization,  and, the  changes  in  the  mean 
height  of  the  ionosphere,  the  estimated  accuracy  of  TEC  data  obtained  from  low  orbit 
satellites  is  not  better  than  +10  percent.  Due  to  the  limited  time  sampling  of  TEC 
from  low  orbit  satellites,  only  seasonal  average  data  can  usually  be  obtained. 

Continuous  temporal  coverage  is  possible  with  TEC  data  derived  from  VHF 
beacon  signals  transmitted  from  geostationary  satellites.  No  loss  of  data  occurs 
during  periods  of  severe  amplitude  scintillation  if  a  good  polarization  follower 
(polarimeter)  is  used.  However,  with  the  Faraday  method  of  TEC  data  determina¬ 
tion,  the  TEC  is  measured  only  up  to  a  vertical  height  of  approximately  2000  to 
3000  km,  depending  on  the  density  at  those  heights.  The  TEC  between  2000  to 
3000  km  and  the  satellite  height  of  36, 000  km  can  equal  the  density  below  this  height, 
at  least  at  nighttime,  and  can  be  10  to  20  percent  of  the  complete  total  in  the  daytime. 
The’ estimated  accuracy  of  determining  TEC  in  the  lower  2000  to  3000  km  by  the 
Faraday  technique,  from  VHFsignals  transmitted  from  geostationary  satellites,  is 
plus  and  minus  5  percent. 


9-2.  TIDt  -  A  LIMITING  FACTOR  IN  OBTAINING  TRUE  TEC  VALUE 

Even  if  a  true  group  <path  delay  measurement  of  negligible  error  was  made, 
along  one  direction  from  a  particular  station,  a  TEC  model  will  still  be  required  to 
correct  that  true  TEC  value  for  use  along  another  slant  direction  from  the  same 
station.  It  is  likely  that  such  a  model  would  be  able  to  correct  for  the  normal 
temporal  and  geographic  variations  in  TEC.  However,  the  limiting  factor  is  the 
existence  of  small  scale  travelling  ionospheric  disturbances  (TIDs)  which  can, be  of 
a  size  smaller  than  the  distance  in  the  ionosphere  between  the  two  slant  directions 
viewed  from  a  single  station.  The  occurrence,  size  and  location  of  a  specific  small 
scale  TID  cannot  be  predicted.  Hence,  these  TIDs,  which  are  typically  from  1  to 
5  percent  of  the  TEC,  may  become  the  ultimate  limitation  to  TEC  corrections. 


63 


Ionospheric  models  which  are  required  to  independently  correct  for  time  delay 
with  a  resultant  50  percent, rms  error  are  relatively  simple  and  are  at  hand.  Three 
of  the  TEC  models  outlined  in  this  report  have  been  tested  against  TEC  data,  and, 
in  the  mid-latitudes,  have  obtained  an  approximate  75  percent  rms  ionospheric 
correction.  However,  a  system  which  requires  an  independent  correction  for 
99  percent  (rms)  of  the  ionospheric  time  delay  had  better  not  be  built!  A  50  percent 
rms  correction  is  easy  to  accomplish  in  most  cases;  99  percent  is  impossible;  the 
tasktof  ionospheric  workers  in  this  field  is  to  get  as  close  to  99  percent  as  is 
practical,  in  the  most  efficient  way. 


9-3.  IMPROVING  PREDICTION  ACCURACY 

The  work  which  should  greatly  improve  prediction  accuracy  involves  the  times 
when  the  greatest  differences  from  median  conditions  occur.  These  times  are 
usually,  but  not  always,  during  magnetically  disturbed  conditions.  An  understanding 
of  the  physical  mechanisms  which  produce  the  observed  changes  in  TEC  during 
magnetic  storms  is  a  long  sought  after  goal;  but,  unfortunately,  there  is  not  yet 
even  a  completely  clear  experimental  picture  of  the  changes  which  occur  in  TEC  on 
a  statistical  basis  during  magnetic  storms  over  all  regions  of  the  earth.  As  more 
data  are  being  collected  a  clearer  picture  will  emerge  and  there  is  reason  to 
believe  that  within  a  very  few  years  a. much  better  understanding  of  the  world-wide 
behavior  of  total  electron  content  will  emerge. 


